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SUMMARY: The programming of the thrust direction of a rocket missile, to 

obtain maximum range over a flat earth by the most favourable employment 

of the aerodynamic forces, is considered for the special case when these forces 

are small by comparison with the missile weight. The cases of a V.2-type 

missile, (a) with lift and (b) without lift, are dealt with in detail and numerical 
results have been computed. 


1. Introduction 


In an earlier paper”, the problem of achieving maximum range with a rocket 
missile by appropriate choice of the thrust direction programme was solved, in 
principle, under very general conditions. In the same paper, a complete solution 
was found for the special case of a missile being fired over a flat earth and subject 
to no air resistance. It was found that the direction of thrust had then to be kept 
fixed. The effect of earth curvature and the consequent variation in gravity have 
been discussed in a second paper. In the present paper, an attempt is made to 
assess the effect which the presence of an atmosphere will have upon the optimum 
trajectory. In general, the calculation of the trajectory of maximum range has been 
shown” to involve the estimation of the conditions at “all-burnt” and an integration 
backwards along the missile track to the instant of firing. The estimated conditions 
at “all-burnt” have then to be adjusted until the rocket co-ordinates and velocity 
components at the launcher are correct. Having programmed an electronic 
computer for the backwards integration, this process of adjustment could be carried 
out without too great difficulty.. However, the author’s resources do not extend 
further than a desk machine and accordingly he has been forced to restrict himself 
to a special case for which it is possible to avoid the laborious procedure, but 
which will, none the less, reveal the main features of the atmospheric effect. This 
is the case for which the aerodynamic forces are small by comparison with the 
missile weight, so that the atmospheric effect can be regarded as a perturbation 
superposed on the known optimum trajectory in the absence of aerodynamic forces. 
Aerodynamic forces would be relatively small for a large intercontinental missile or 
for a missile launched at high altitude, though possibly not quite small enough to 
Satisfy fully the assumptions of the present analysis. However, it may be expected 
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that the main characteristics of the optimum thrust direction programme, in the 
case where the aerodynamic forces are relatively large, will be qualitatively similar 
to those for the special case to be considered, and hence that the results obtained 
will be of significance under more general conditions than those assumed in their 


calculation. 


NOTATION. 


see equation (7) 


values of 2,,, 8, respectively on the optimum trajectory in the 
absence of aerodynamic forces (this will be referred to as 
trajectory A) 


first-order corrections to A,,., B,, to allow for aerodynamic forces 
drag component 
acceleration component due to motor thrust 


components of missile retardation due to gravity and/or aero- 
dynamic forces 


values taken by the expressions for F,, F., on trajectory A 
acceleration due to gravity 

see equation (4) 

direction cosines of direction of motor thrust 

lift component 

values of /,. /, on trajectory A 

first-order corrections to L,, L, to allow for aerodynamic forces 
rocket mass 

see equation (27) 

see equation (28) 

rocket’s co-ordinates at “all-burnt” 

rocket’s velocity components at “all-burnt” 

values of q,, q, on trajectory A 

first-order corrections to Q,, Q, to allow for aerodynamic forces 
values of g,, G2 on trajectory A 

first-order corrections to O.. OQ, to allow for aerodynamic forces 
missile range 

time after launching 

value of ¢ at “all-burnt” 

rocket velocity 


co-ordinates of the rocket 
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X,,X, values of x,, x, on trajectory A 
2X,’,2X,’ first-order corrections to X,, X, to allow for aerodynamic forces 
Yin, Zix See equations (5) 
Yi, Zx Values of On trajectory A 
eY,,’,¢Z;, first-order corrections to Yj, Z;, to allow for aeredynamic forces 
By, see after equation (7) 
8 angle of yaw 
8 Kronecker deltas 


m 


small parameter multiplying aerodynamic forces 

6 inclination of the rocket’s velocity to the horizontal 
® inclination of the rocket thrust to the horizontal 

p atmospheric density 


pf) atmospheric density at sea level 


2. An Extension of the General Theory 


If Ox,, Ox. are horizontal and vertical axes respectively through the launcher 
and lying in the plane of the trajectory, the equations of motion of the missile are 


(i=1,2) . ‘ ‘ (1) 


where f is the acceleration due to the motor thrust, (/,, /.) are the direction cosines 
of the direction of thrust and (-F,, - F.) are the net components of force per unit 
mass acting upon the rocket due to all agencies but its own thrust. In Ref. | it was 
assumed that the F; were explicitly expressible as functions of the variables x;, x;,M 
(rocket mass) and t. However, since the aerodynamic lift force will depend upon 
the orientation of the rocket axis relative to the direction of motion, the theory will 
now be extended by assuming that the direction of this axis is along the line of 
thrust and hence that 


Proceeding as in Ref. 1, it can now be shown that the thrust direction 
programme for maximum range R is determined by the equation 


OR OR 


where Qu+ =: . ; (4 
q a4, 


*The summation convention is operative in equation (4) and throughout the papet. 
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and (q,. q2), (G,.42) are the rocket co-ordinates and velocity components respec- 
tively at “all-burnt.” The functions z,,, 8,, have already been introduced in Ref. 1 
as follows. 


Let G=1,2;. £=1,2, 3, 4) ‘ (5) 


represent a linearly independent set of solutions of the equations 
resp 
ex, xj . . & itw 
Forming the determinant 
Yun Vis Yas Vue | 
Yor Yoo Yo3 
the quantities z,,, 8,, are the determinants formed from A by replacement of its 
s™ row by the elements 
Yri (T), Yre (T), Yrs (T), Yrs (T), 
Sub 
3. Solution for Small Aerodynamic Forces 
Neglecting earth curvature and rotation, the equations of motion (1) will be 
written in the form 
x, fl;, (8) whi 
Bt 
where (-:F,, -¢F,) are net components of the aerodynamic forces per unit mass 
and < is a small parameter. whi 
Sati 
If, now, in the absence of an atmosphere, the optimum trajectory and thrust 
direction programme are determined by the equations 
at 
May 
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respectively, the presence of the small aerodynamic forces can be allowed for by 
replacement of these equations by 


x= X, (+ 2X/ () +O 
L=L,(t)+ (t)+ O (¢), ? (12) 


respectively. Substituting from equations (11) and (12) into equations (8), and 
equating coefficients of like powers of = from the two members of each equation, 
it will be found that 


X/+F/=fL/ (i=1.2), 14) 
where F;’ is the result of substituting x;—X;, /;=L; into F;. 


Similarly, it will be assumed that 


(15) 
+ 2Z/ +O 
Substitution in equations (6) then yields the equations 
Z<-0, . . . . @ 
OF/ 
Zi + —- Y;=0, 
ax, (17) 
Y/-Z/=0 


If the fundamental set of solutions of equations (6) is taken to be the one 
which at t=T (“all-burnt”) satisfies the conditions 


Zin =Sins Vir = bs (18) 


where 5,,=0 (i ¢ k) and is unity otherwise, then, since the conditions (18) are to be 
Satisfied for all «, 


at t=T. 
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The four sets of solutions of the equations (16) satisfying the conditions 


(19) are pl 
Gh Z,=0. Z,,=6. Y,,=1. Y,,=0; (21) 
(iii) Z,,=0, Z.3=0, Y,;=1, Y,, =0; 
(iv) Z,,=0, Z.,=0. Y,,=0, Y.,=1 
Employing each of these sets of values for the quantities Y; and Z; in equations N 
(17) in turn, corresponding sets of values for the quantities Y;’, Z,/ may be found by 
integration backwards from t=T, subject to the initial conditions (20). The 
quantities 0F;’/0X,, OF ;’/0X; are to be calculated along the optimum trajectory in 
the absence of an atmosphere. This trajectory can be found by the method of 
Ref. 1 and thus these quantities are available for integration. 1 
The determinants z,,, 8,, may now be computed in the form a 
il 
2A,, +O (?), 
(22) 
+ +O (€*), 
where A,,=T-t, A,,=90, 
A,,=0, A,.=T —t, 
(23) 
B,,=1, B,,=90, 
=0, = 1, 
A,,/= Yo. + (T -1) (Y,,'+2Z,,’) (T - (24) 
} 
B,./=Y 
At “all-burnt,” suppose that 
=Q;+2Q/ +O 
(25) 
Gi=Qi+ +O (2), 
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ons Then, from equation (3), it will be found that the optimum thrust direction 


programme is determined by the equation 


, 


2 
i, L.+2L,’+O(8) (26) 
1) 
where Ni= - on - By. (27) 
a0, ad, 
OR 4 , OR 
ns N/= - + A; i R 
20,0. 20,00. 
The quantities Q,, 0; specify the missile position and velocity at “all-burnt” in the 
absence of an atmosphere and can therefore be found without difficulty. Accord- 
ingly, the quantities N;’ can be computed as linear forms in the unknowns 
Q.’, QO.’ from the equations (28). 
2) Multiplying across in equation (26) and equating coefficients of like powers 
of « from the two members of the resulting equation, it will be found that 
L,N.=L.N,. 
But /,7+/,?=1 identically, and hence 
Equations (29) and (31) now imply that 
+) 
and equations (30) and (32) that 
NAN.N,’ N,N,’) N(N,N,’ - N.N,’) 34 
Equations (34) determine the L,’ as linear expressions in the quantities Q,’, ay. 
Substituting these expressions in the right hand members of equations (14), the 
) & quantities X/ are also given as linear expressions in these unknowns. Integration 
of X/ with respect to f over the interval (0, 7) yields QO’. and a second integration 
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over the same interval yields Q,. Thus, the quantities Q/, Q/ can be expressed 
linearly in terms of themselves, thus providing four simultaneous equations from 
which their values can be computed. 


The N; are now completely determined and the thrust direction programme 
can be found to O (:). 


4. Application of the Theory to a V.2-type Missile 


The theory of the previous sections will now be applied to the case of a missile 
having aerodynamic characteristics approximating to those of the V.2 missile. Since, 
for such a missile, the aerodynamic forces are not small by comparison with the 
weight. it will be supposed that they are reduced by a factor : by, for example, 
supposing the atmospheric density reduced by the same factor. An equivalent 
assumption would be that the thrust and weight forces are increased by a factor 1/<. 
The effect on the thrust due to the difference between jet and atmospheric pressures 
at the venturi exit will be neglected. Units will be the usual absolute units of the 
f.p.s. system (i.e. forces will be expressed in poundals) unless otherwise stated. 


It will be assumed that the aerodynamic forces acting upon the missile can be 
resolved into a drag component D acting in a direction opposite to that of the 
motion and a lift component L acting in the upwards sense and perpendicular to D. 
Then, if # is the angle made by the direction of motion with the x,-axis, 


F,=(Dcos 6+L sin @)/M, 


F,=(Dsin6-Lcos@)/M, 
and tan 6=%./x,. . (36) 

If « is the rocket velocity, so that 


p is the atmospheric density and p, is the density at sea level, the drag will be 
supposed given by 


Assuming a steady temperature lapse rate of 2°C per 1,000 ft. and taking the 
launcher to be at sea level, then 


0 
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If » is the angle made by the rocket axis with the horizontal, so that 


cos @=1,, . ‘ . (40) 


then = —46 (radians) is the angle of yaw. The lift will be supposed determined 
by the formula 


L=6-28 


Taking the initial mass of the missile as 28,500 Ib. and the rate of propellant 
expenditure as 275 Ib./sec., the rocket mass at time ¢ is given by 


M = 28,500 275t. (42) 

Also, if the motor thrust is taken to be 60,000 Ib. wt., then, with g= 32:2, 
f =1-932 x 10°/M. (43) 
Equations (35), together with equations (36)-(41), now determine the F; as 


functions of the quantities x;, x;, |, M, t. The partial derivatives of the F; with 
respect to the first six of these variables are: — 


OF, OF, _ 
Ox, Ox, 
= fos cos 26) D+ L sin 26+ cos 28) | /(2Mz), 
[p sin 26 L (3.+cos 28) - sin 26 | /(2M~), 
= [p sin 20 + L (3 cos 26) - sin 26 | /(2M»z), 
- (44) 
[p (3 cos 26) - L sin 26+ (1. +cos 2#) | /(2Mv), 
-0-3918 x 10-* 4, 
y 0-3382 x 10 5 sin 6, 
1 OF, 
f al, =0-3918 x 10 5 008 4, 
1 OF, 
= ~0°3382 x 10-* ‘cos 
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A thrust duration of T=70 sec. was assumed. After the cessation of motor 
thrust, the effect of the aerodynamic forces upon the subsequent motion of the 
missile was neglected. Thus the horizontal range of the missile from “all-burnt” 
can be calculated from the in vacuo formula 


(4.4 +2gq,)}. 


The net range is then given by 


The optimum trajectory for the missile in the absence of aerodynamic forces 
has been calculated in Ref. 1. In this case, the direction of thrust is constant and 


The values of the partial derivatives (44) at points along this trajectory were 
computed (five-second intervals) and, by a process of numerical integration, the 
values of the quantities Y;,’, Z;,’ were found at these same points from equations 
(17) and initial conditions (20). 


Substitution in equations (24) now provided tables of values of the A;;’, B,,’. 


The conditions at “all-burnt” on the optimum trajectory in the absence of 
aerodynamic forces were taken from Ref. 1 to be as follows: — 


Q,=147,700, Q,=92,250, 
QO,=5,171. QO, = 3,737. 


The numerical values of the partial derivatives 0R/0Q,, 0*R/0Q;0Q);, and so on, 
were then computed and the quantities N;, N; tabulated from the equations (27) 
and (28). Tables of values of the L,’ now followed from equations (34) and a 
further four numerical integrations of the equations (14) resulted in simultaneous 
equations for the unknowns Q/, Q/. as explained in the previous section. The 
solution of these equations proved to be 


(48) 


Q,’ =882,000, Q.’ = - 618,000, 
Q,’= 6,850, Q.’= - 930. 

The functions X,’ were now tabulated, thus determining, to the first order in <, 
the effect on the optimum trajectory of the aerodynamic forces. The optimum 
trajectory in the absence of these forces is shown, drawn to scale, in Fig. 1 as the 
curve A. Taking «=0-05, the optimum trajectory in the presence of these forces 
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is shown, to the same scale, as the curve B. With such a large value for ¢, it is 
probable that second-order effects are significant, but this value has been taken with 
the object of magnifying the divergence from the in vacuo trajectory, in order that 
its nature will be clear from the figure. The direction of thrust (and hence of the 
rocket axis) is indicated by arrows at 10-second intervals. 


It will be observed that the early portion of the track B is approximately 
horizontal (in fact, it dips below the earth’s surface initially, although, of course, 
this cannot correspond to a practicable trajectory). During this time, the missile 
gathers momentum and its angle of yaw increases to provide a steadily increasing 
lift force, by means of which it eventually gains altitude. During the latter stages 
of the motion, the angle of yaw. is decreased. Thus, broadly speaking, the lift 
assists in the generation of the vertical momentum permitting the rocket thrust to 
be employed primarily in the generation of the horizontal momentum. Under 
normal atmospheric conditions it may be expected that, for the optimum trajectory, 
the lift will be responsible for the greater part of the vertical momentum. 


The manner of variation of the angles » and 4 with the time is shown 
graphically in Fig. 2. 
The increase in range, beyond that attained in the absence of the aerodynamic 
forces is given by 
oR OR 
Q/+ . OF ) =1-64 x 10°. 49) 
30.2 20, Q ) ( 
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It is interesting to observe that the range increment is positive and, with <=0-05, 
amounts to about 16 miles. If, however, the action of air resistance during the 
ballistic phase is taken into account, the apparent increment may well be converted 
into a decrement. Nevertheless, it is clear that, by appropriate choice of the thrust 
direction programme, the adverse effect of drag upon the range can be offset, to a 
large extent, by a compensating effect provided by the lift. 


As a check upon the validity of these conclusions, the whole problem was 
recomputed, neglecting the lift force. The optimum trajectory for a V.2-type 
missile subject to drag alone was therefore found. The following values were 
obtained for the quantities Q,’, Q/, 


Q,’= 104,700, = 53,800, (50) 
O,’= 2,050, O,’—240. 
With these values, it will be found that 
dbR=—4-9< x 10°, . (51) 


i.e. there is a range decrement which amounts to about 5 miles when «=0-05. 
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The optimum trajectory for this missile (with :=0-5) is shown as curve C in 
Fig. 1. In this case, no advantage is to be derived by remaining in the lower strata 
of the atmosphere. Thus, the thrust is now employed to gain altitude fairly rapidly, 
in order that advantage may be taken of the lower atmospheric density, and 
consequent lower drag, at the greater heights during the phase of acceleration 
towards maximum velocity. The thrust direction is again shown by arrows at 
10-second intervals. 
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The Compressible Flow Through Cascade 
Actuator Discs 


J. H. HORLOCK 
(University Engineering Laboratory, Cambridge) 


Summary: A theory of the incompressible flow through two- and three- 
dimensional cascade actuator discs has been developed by several workers over 
the past ten years, and its accuracy has been confirmed in several experiments. 
This theory is briefly reviewed, and a parallel theory for subsonic compressible 
flow through actuator discs is developed. Approximate solutions for several 
examples are considered, including a compressible shear flow through a two- 
dimensional cascade, and a compressible flow through an annular cascade 
of guide vanes. 


1. A Review of Incompressible Flow Through Cascade Actuator Discs 


1.1. INTRODUCTION 


The incompressible flow through cascade actuator discs has been studied by 
several authors”. A brief review of approximate solutions to the two- and 
three-dimensional problems is given here, to provide a basis for comparison with 
the compressible solutions derived later. 


Referring to Fig. 1, a co-ordinate system x, v, z is considered, with velocities 
u,v, w and vorticities €, 1, ¢. 


CO-ORDINATE SYSTEM 
AND NOTATION 


FIGURE 1. The  two- 


dimensional problem. ACTUATOR— 
PLANE 


STATIONS 02 2 
Received July 1957. 
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CO-ORDINATE SYSTEM 


AND NOTATION wit 
ACTUATOR —— 
PLANE 
FiGURE 2. The _ three- 
\ dimensional problem. 
x20 x 
STATIONS ' ol 02 2 


A cascade actuator disc made up of an infinite number of blades each of very 
small axial length, is located at x=0. The flow from a station (1), far upstream, to 
a station (2), far downstream, is considered. At both (1) and (2) the static pressures 
will be constant. The cascade acuator disc will provide a discontinuity in tangential 
velocity v and vorticity », but there is no change in axial velocity and radial velocity 
across the disc. Since there are many blades, all derivatives in the y-direction 
vanish. 


Ruden™ originally made the two assumptions that the tangential velocity and 
vorticity changed little in the x-direction (i.e. v, »=f (z)) for a “small perturbation” 
flow. Hawthorne and Armstrong (for two-dimensional cascades) have shown that 
these two assumptions are valid if w/u, 1, 0n/0z, dv/0z are always small. In three- 
dimensional flow (see Fig. 2), these conditions become that w/u, », 0n/dr, 0 (rv)/é@r 
should be small. 


Further, if the tangential vorticity is constant along lines of constant z or r, then 
a consideration of the axial velocity fields induced by the 7-component of vorticity 
leads to the conclusion that in incompressible flow the axial velocity at the disc is 
the mean of the axial velocities (at the same z or r) far upstream and far down- 
stream. In the notation of Figs. 1 and 2, 


Uy, = Uo. = (u, + 


These assumptions were found to be reasonable in an exact analysis for two 
particular examples considered by Bragg and Hawthorne”. 


NOTATION 


x,y,Z co-ordinates 
u,v,w velocities in x, y, Z co-ordinate system two-dimensional flow 
£,,¢ vorticities in x, y, z co-ordinate system 
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| length of actuator disc in z-direction 
A., constants 
x,9,r co-ordinates 
u,v,w velocities in x, 4, r co-ordinate system three-dimensional flow 
vorticities in x, r co-ordinate system 
stagnation enthalpy 
a air angle 
p density 
p static pressure 
a_ speed of sound 
M.. axial component of Mach number 
A= /(1-M,,.”) 
area 
X=4 (Uy Usm) (U, — 


M, tangential component of Mach number 


1- Ma." 

1+Mo,. 
1 


q= V (u? 
W flow rate 


Subscripts 
1 axial location far upstream of an actuator disc 
2 axial location far downstream of an actuator disc 
Ol axial location immediately upstream of an actuator disc 
02 axial location immediately downstream of an actuator disc 
x axial location far upstream or downstream of an actuator disc 
t tip 
h hub 
m mean 
Superscripts 
perturbation 


* relating to conditions where M=1 (in Tables IV, V and VI). 
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1.2. ANALYSIS FOR Two-DIMENSIONAL INCOMPRESSIBLE FLOw 


The continuity equation for two-dimensional incompressible flow is 


Ou ov Ow _ 


and, as derivatives with respect to y must vanish, 


du 


ox az (1a) 


If the velocity u is made up of the velocity u,, (z) at the location far away from 
the disc and a small perturbation w’ (u=u,,+u’), then equation (1a) becomes 


Ou 
- + 


Ox 0, 


where wu’ and w are both small compared with u,.. 


The Ruden assumption for the vorticity, 


N= Noo (2), 
Ou,, , Ou’ OW 
may be written — eae 
ou’ = Ow 


From equations (1) and (2), 


= 
a solution of which is 
w= Se ***[A, sin kz + B, cos k,z). (4) 
i 


The boundary conditions for this equation are that w should vanish : — 
(i) at the side walls (z=0, J). 
(ii) far upstream of the disc (for the flow upstream). 


(iii) far downstream of the disc (for the flow downstream). 
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With these restrictions 


w= = Aww sin (nzz/De"*'' upstream, 
(5) 
w= Ang Sin downstream. | 


Integrating (1a) with respect to x, from infinity to x, gives the solution for vu’, 


u’= - Ay, cos upstream, (6) 
u’= + Ang cos downstream. | 


The second Ruden assumption, that u,,=u,.=4 (u, +u,), provides a boundary 
condition for wu’ at the disc 


Uy,’ =4 — U,)= — & Anu COS | 


Uys’ = -4(u,—u,)= + Ana COS (n%Z/1). 

If the distribution of uv, and u, with z is known, then the coefficients A, may 

be obtained from a Fourier analysis. For most cases met in turbo-machines, 

since 4(u,—u,) will be approximately zero at z=//2, and approximately as much 

positive at one wall as negative at the other, the first Fourier coefficient is much 

larger than the later terms. For a linear variation of 4 (u, — u,) with z, for example, 

the second term is one ninth of the first, and decays at the rate e~>**/' away from 

the disc. Thus a solution using the first term only will be sufficient for most 
practical purposes and 


uw =~4(u,—u,)e*'! upstream, | (8) 
u’ = —4(u,-u,)e~*' downstream. 


This solution is exact if 4 (u, - u,)= A, cos (=z/I). 


1.3. ANALYSIS FOR THREE-DIMENSIONAL INCOMPRESSIBLE FLOW 


In three-dimensional flow, equations corresponding to equations (1b), (2), (3). 
(4), are 


9 
ou’ ow 
1 ow w 
w= Se**"[AJ, (kix) + BY, (kv)]. Phase 
i 
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The boundary conditions that w=0 at r=ry, r: give 
J, (kite) (krd=J, (kind Y, (Kirn), 


Table I gives solutions obtained from Jahnke and Emde“*. 


TABLE 
k (r,—ry) 
1-0 
15 3-1609 6-293 


25 3°1966 6°312 


Taking only the first term in the solution again, and using the Ruden 
assumption for the velocity at the disc 


u’ =4(u,—u,)e™ upstream, 
2 3 1 | ( 14) 
u’=—4(u,-u,)e “ downstream, 


where k, is given in Table I, as a function of r,/ry. 


1.4. THE DIRECT AND INDIRECT PROBLEMS IN INCOMPRESSIBLE FLOW 


So far only the indirect problem of the flow in turbo-machines has been con- 
sidered—the flow in a field in which u,, u. are specified far upstream and far down- 
stream. In the direct problem the flow far upstream, and a distribution of outlet 
angle 2,. at exit from the disc, may be specified (where tan z).=7,./U 2). It is 
required to find the velocity distribution far downstream and elsewhere in the field. 


A solution to this problem is obtained by using the equation of motion in the 
r-direction (for three dimensions) in the form 


where H is the stagnation enthalpy and 


For a stationary guide vane row, for example, accepting a uniform axial velocity 
profile far upstream, H,,=H,,. 


2H ys _ 

or 
Vor O 
r Or 
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Since 

or 
Ou, + Vor 0 
Or = 


But = Ugo tan Aye 


Ou, 


and 


K 


du, 


ar 


(rv,.)=0. 


ap (ru,, tan 2,.)=0. 


For the particular case in which tan z,,=constant, using the Ruden expression 


for Uy., 
dus (4 tan’ v2) ls tan*2,,, 


It is of interest to note that the solution to 


this problem based on a “radial 


equilibrium” approach”, by assuming there is no radial velocity at exit from the 
blades (i.c. that station (2) is at the trailing edge of the blade row where the angle 


distribution z., is constant), would give 


on, @ 
or r or TH, 


and ur*'™. — constant. 


ou 
or 


2 =0, 


(17) 


Considering the flow through a set of guide vanes where the exit angle is 45° 


and r,/r,=0°4, equations (16) and (17) become 


or 


constant, 


u.r'/* =constant. 


(16a) 


(17a) 


Table Il shows a comparison of root and tip axial velocities obtained from 


these two solutions. 


TABLE Il 
from (16a) 
r,=0°4 100 120-5 
n=10 100 88°5 
116 


= 


from (17a) from (16a) 
131-8 110-2 
83-2 94:2 
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To illustrate the indirect problem of specifying an angle distribution to 
produce a particular velocity distribution, the design of a guide vane row to produce 
a constant tangential velocity distribution 7, may be considered. 


Far downstream the axial velocity distribution is obtained from equation (15), 


du, v.d 
“sa (rv,)=0. 


and u.= V(A-2v.?logr), (18) 


where A and v, are constants. 


The angle distribution at a cascade actuator plane where the axial velocity is 
(u,+u,) would be 


y 2v? 


tan ¢,.= 


but the angle distribution at the trailing edge of a blade where the flow is in radial 
equilibrium would be 


tan z,— (20) 


V(A -2v.? log r)’ 


A comparison for a guide vane row producing a tangential velocity 7, equal to 
the axial velocity u,, far upstream, is given in Table III. 


TABLE Ill 
uy | | v, | | Bos 
| | | from (20) from (19) 
r,=0°4 100 146 100 34-4 39-1 


r,=10 100 53 100 62:05 52°55 


1.5. EXPERIMENTAL WORK IN INCOMPRESSIBLE FLOW 


Several experiments, on two-dimensional cascades and axial flow compressors 
have been made“: '”, all at low speeds. These results confirm the accuracy of the 
analysis given in the previous section. For example, Figs. 3 and 4 show axial 
velocity profiles observed far upstream, at the leading edge, at the trailing edge, and 
far downstream of a single stator (Fig. 3) and a single rotor (Fig. 4). The exponential 
decay in axial velocity up to and away from the blade row is clearly seen in each 
case, and the velocities are compared with calculations made from the analysis 
outlined in the previous section, using cascade data to give the outlet angle distribu- 
tion. In each case, for the purposes of calculation, the actuator disc was supposed 
to be located on the axial centre line of the blade row. 
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Axial velocity profiles. 
FiGurE 3. Stationary row. Ficure 4. Rotating row. 


2. Compressible Flow 
2.1. ANALYSIS FOR Two-DIMENSIONAL COMPRESSIBLE FLOW 


A similar approach may be made to the problem of compressible flow through 
the two-dimensional system of Fig. 1. 


The continuity equation is 
0 
ax (pu) + By (pv) + (pw)=0, 
but, since the y derivatives vanish, 


0 


The Euler equations of motion are 


lop_ ou du Ou 

1 dp ow Ow ow 
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Noting that the y derivatives again vanish, and that the velocity of sound 
a= ¥ {(dp/dp),}, these equations may be rewritten 


a’ Op Ou eu 


ov v 
0= + (23a) 


a op Ow Ow 


Substituting in equation (215) from equations (22a) and (24a), 


(1- + ow) = =0 . . (25) 


If small perturbations from the flow at +0 are considered, then u=u, (z)+u’, and 
u’/u,, wW/u,, are both small. Further, if another of the assumptions of the incom- 
pressible analysis is made, that 7, is always small, then terms such as 0u/0z, Ow/0x 
are small, and equation (25) becomes 


ou’ ( ‘ ow 
2 = 2 
a 1 M.. ) + az 0, (26) 


where M,=u/a is the axial component of the Mach number of the flow. 


The variation of v with x is seen to be of second order, from the second Euler 
equation, for dv/d0x= -—(w/u) dv/0z, and both w/u and 0v/0z are small. 


The variation of » with x may also be shown to be of second order. The second 
Helmholtz equation may be written 


dn 


ou 


But €= dv/0z and (=dv/0x, and equation (27) becomes 


on won dw 7 
dx u dz az + 


Since both » and @1/0z are small, the variation of » with x may be neglected. 
Using this result, 


Ou,, Ou Ow 
+ 


Ou Ow 
29 
(29) 
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From equations (26) and (29), the equation 


ae (1 dz + az =0 (30) 


is obtained. 
Writing the axial Mach number 


equation (30) becomes 


Since = 0u,./0z is small, du’/0x and 0M,./¢z will be small, and equation (30a) may 
be written 


ax! (1 - + dz? 0. . (30d) 


A solution of this equation is 


w—- Seth {A, sin {k,/ (1 — M,.*)} z+ B, cos (1 - M,..*)}z]. 
1 


With the boundary conditions w=0 at z=0, /, 


where A= /(1 - M....”), and the solution becomes 


w= sinnzz/! upstream, 
1 


(31) 
w= Ye sin (nzz/1) downstream. 
1 
ou’ , 
Equation (26), which is + =0, 
may be integrated from the infinity station to x and 
w= — (4,,/A,) cos upstream. | 
(32) 


wW=+ Le cos downstream. 
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2.2. THE AXIAL VELOCITY AT THE DISC IN COMPRESSIBLE FLOW 


To find the boundary conditions on axial velocity at the disc, the Ruden 
assumption must be reconsidered for compressible flow. 


The induced velocity due to a vortex field, », (z)=du,/dz, stretching from 
x= - tox= +00 is at any point x (including x=0) where u,,, is the mean 
axial velocity across the section from z=0 to /, i.e. u,.=( ‘u,d2)/I. The induced 
velocity at x=0 due to a vorticity field »,(z)=du,/dz, stretching from x= - © to 
x=0 will be 4 (u, 


Similarly the induced velocity at x=0 due to a vortex field 7, (z)=du,/dz, 
stretching from x=0 to x= +00 will be 4 (u.- Won). 


The total induced velocity at the disc (x=0) due to the combined fields 
(z)=du,/dz (- 0 <x <0) and ».(z)=du,/dz (0< x < will be 4(u,+u.)- 
$(Uym+Uom). 


Thus the axial velocity at station (01) is 
Uy, = Uym +4 (Ue + U,) — + (Un +,)+4 (ym om), 


and at station (02) 


The change in axial velocity across the disc is 


Mog — Ug, =Uom — Uym. 


and is constant with z. 


In incompressible flow u,,=Uy2. SINCE U4 =Uym, but in compressible flow u,,, will 
be greater or less than uv. according to whether the row is a compressor or a turbine 
type row. The perturbation velocity at station (02) is 


=4 =x... . (34) 
These boundary values may be used with the first term only of equation (32) to give 


w=xe upstream, | 
(35) 


u’=-xe downstream. 
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2.3. DENSITY VARIATION IN COMPRESSIBLE FLOW 


Writing the density p=p.(z)+p’, and ignoring the cross products of small 
terms, equation (21b) becomes 


ap’ Ou’ aw 


_ 
ax + Po ax + +Ww =(), (36) 


u Pao dz 

For two-dimensional cascades, the pressure at the infinity station is constant. 
and dp,./0z will be zero if the static temperature is constant. However, even if 
small variations in the static temperature and density p,,(z) are allowed, the last 
term of equation (36) is of second order of smallness, and the equation becomes 


ax ax <7 


Substituting for 0u’/dx and dw/0z from equations (31) and (35), and integrating 
from the infinity station to x, 


(37) 
+ Pow M.?e-** downstream. 
2 
(1 - @kx ) upstream, 
(37a) 


M.” 
p=p.(1 + ts) downstream, 


2.4. EXAMPLES OF Two-DIMENSIONAL COMPRESSIBLE FLOW 


As one example of the compressible flow through a two-dimensional cascade, 
a uniform entry flow may be considered (u,, v, =constant). If the cascade exit angle 
is also uniform, the exit velocities u,, and u, will be uniform. There will be no 
perturbation in axial velocity or density upstream or downstream of the disc, since 
X=4 Solutions for u,, and u, based on actuator disc 
theory or “radial equilibrium” will be the same. For example if M,=0-5, 2,=0, 
2, = 30°, a solution may be obtained directly from gas tables“. Table IV (p. 126) 
gives such a solution, using the notation of Ref. 11. 


As a second example of a two-dimensional compressible cascade flow, a linear 
variation of velocity and Mach number at entry to a cascade (a,=0, %2=30°) 
is considered. 
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Writing the equations of motion in the z-direction, 


6) Wes du, 
Voo az + 


== Me, ‘ 
st 
since 
= 
it follows that 


2 (p. tan?’z,,. 
(2+ tan?z,.) 


For the particular example, 


p./p, 0-936. u.=0-65u,+constant, . . (38) 


using the solution of Table IV (p. 126) to give the density ratio p.//,. 


A trial and error solution gives the velocity distribution in the complete field, 
as in Table V (p. 126), in which the constant of equation (38) is taken as 41-1. 


To check the analysis, values of u,, as calculated directly from u,,=u.-— x and 
as calculated from the density variation p,./p, and from Ref. 11 are compared in 
the last two columns. 


2.5. ANALYSIS FOR THREE-DIMENSIONAL-COMPRESSIBLE FLOW 


Equations corresponding to (215), (22a), (23a) and (24a) in three-dimensional 
flow, in which the 6 derivatives vanish, are 


ax (pur) + (prw)=0, . (9) 
_@ op _ 
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Substituting in equation (39) from equations (40) and (42), the equation 


Ou M,2)+ On (1 - +M,2) | 
Ox or r 


ow 
3 
a \ér x) 


is obtained. 


As in the two-dimensional case, if » and 0y/0r are small then »=7,, (r) and 


Ou daw 
44 
or Ox 
Also, since » is small, equation (43) becomes 
Writing M.=M.,, (r)+M,’ and Ms= (r) + Me’. 
and assuming 0M,,./Or, 0M¢,,/0r are small, the equation 
w 
may be obtained by cross-differentiation of equations (44) and (435), where 
Max" 
The solution* is w= Z 2+M, 2 . (46) 


The boundary conditions on this equation are w=0 at r=r, and r=n. 


If Ms,,—> 0, this requirement is the same as that derived in the incompressible 
analysis, 


J, (Karn) Y, J, (kr) (Kary) (13) 
and the solutions of Table I are valid. 


*See Ref. 13, p. 146 for the general solution of Bessel type differential equations. 
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If M, —> 1, then the boundary condition requires that 
Fars (Karn) 4 (kar) Y (Kirn). (47) 
Solutions to this equation are 


1-0 
| - ry) 3:17 


However, it should be noted here that the approximate solution is valid for low 
values of My only. The assumption that dp,,/dr is small has been made, and 


or @ Or oer @ r 


Thus if dp,,/dr is small, then Mo,.* must be small. 


It is clear that the tangential Mach number makes but little difference to the axial 
decay factor. Taking the first term only 


kxivQl M cx 2) 


where l=n-Tr and kl=~7=. 


By writing p=p,,(r)+p’ and assuming dp,,/Or is small, the perturbation density p’ 
is obtained, 


2) 


Use 
p’= downstream, 
M 2 kx/ J/ —M,2) 
and ( 1 - ) upstream, 
(37a) 


kx/V(Q M.,2) 
p= ) downstream. 
2 


2.6. EXAMPLE OF THREE-DIMENSIONAL COMPRESSIBLE FLOW 


As an example of the three-dimensional case, the indirect problem of design 
of a guide vane row to produce a constant tangential velocity is again considered, 
but this time in a compressible flow of entry Mach number 0-5. 
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The analysis for the incompressible flow is valid for the velocity u.,, 


u,= /(A logr), . . (8) 


but if v.=100, as before, the constant A must be altered because of the drop in 
density through the guide vanes. Table VI shows an approximate solution; again 
values of u,, obtained from direct calculation and from the density calculation and 
from Ref. 11 are compared. 


3. Discussion 


The approximate solution for compressible flow should be limited to flows in 
which the axial Mach numbers do not approach unity, and the tangential Mach 
numbers are small. However, even in the example considered, of the compressible 
flow through a row of guide vanes, where the maximum tangential Mach number at 
exit is approximately 0-5, the density perturbations given by the analysis appear to 
be consistent with the velocity perturbations predicted. (It is of interest to note 
that the pressure variations may be transmitted axially upstream if the axial Mach 
number is less than unity, even if the absolute Mach number exceeds unity; but as 
the solution given has been derived for axial Mach numbers less than unity and 
small tangential Mach numbers, it would not be expected that the solution would 
be accurate in the example given in Table VI). 


The particular case of irrotational flow through a “free vortex” guide vane row 
is interesting, since it illustrates further the restriction of low tangential Mach 
number and small values of ¢p,,/Or=dp.,/dr. The solution given suggests that 
there will be no perturbations in axial velocity on either side of the actuator disc, 
SINCE UW, =Uym, Ug=Ugm and x= 4 (u, U,) 4 This is understandable 
for there is no tangential vorticity in a “free vortex” flow. Further the density 
perturbation p’= +(xM,,,?/u,,) Pp. must also vanish. Since there is a pressure 
gradient with radius dp,,/dr=p,v,,°/r far downstream, there must be a density 
variation also, and this will be the same as the density variation immediately down- 
stream of the disc. This argument suggests that a radial displacement of the 
streamlines must take place at the-disc, but not upstream and not downstream; since 
the disc is infinitesimally thin, there is a discontinuity in the streamline at the disc. 
However, it should be remembered that it was assumed in the analysis that dp.,/dr 
was small, and any discontinuity in the streamlines at the disc must vanish under 
such restrictions. 


The continuity equation is satisfied at the disc by equations (35) and (37). 
Upstream of the disc, the flow rate W is given by 


r 


t 
W= | dr 


"h 


) rar 
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FiGurE 5. Axial velocity profiles. Example of flow through guide vanes. 


R 


and, since the term involving \* is small compared with the other terms, this 
becomes 


W = Jr, (1 - M.,*) \rdr. 


"h "h 


Since dp,,/dr and dM,.../dr are small; the second integral may be written in 
terms of p,,, and M,,,, the mean values of », and M,, 


Pym (1 Mim?) yr dr = pym - Mim?) rdr- | uot dr rdr+ | dr | 


"h "h 


since, by definition, “= "' 


and similarly W= | p.u.r dr= dr. 
"h 
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FiGurRE 6. Flow through guide vanes. Axial velocity and density changes. 


Finally it is of interest to compare the axial velocity profiles calculated in 
Tables III and VI, for the incompressible and compressible flows through guide 


vanes giving a constant tangential velocity at exit (Fig. 5). In the incompressible 
case there is axial acceleration at the root and axial deceleration at the tip, on both 
sides of the disc. In the compressible flow there is acceleration at the root up to 
the disc, across the disc and downstream of the disc. At the tip there is deceleration 
up to the disc, acceleration across the disc, and deceleration downstream of the disc. 
The compressible flow will result in less radial shift of the streamline, because of the 
increase in pressure and density towards the tip at exit from the row. The changes 
in axial velocity and density with axial distance x// are shown in Fig. 6, for the 


example of Table V. 


4. Conclusions 


An analysis to determine the “axially” subsonic compressible flow through 
cascade actuator discs has been outlined. The decay effects on axial velocity 
observed in incompressible flow should be slightly less in compressible flow, since 
the exponential exponent of such decays is increased by the factor 1/ /(1 - M..,”), 
where M,,,, is the axial Mach number of the flow. Interference effects between 


May 1958 129 


! 
ROOT 
160 
140 
120 
TIP 
TIP 
TIP 
0-7 
o6 ROOT 
+@ 


J. H. HORLOCK 


several rows may be calculated as in the incompressible theory, using the modified 
exponents given, and a design method similar to that of Ref. 12 may be used. 


The density changes up to and away from the disc may be a considerable 


percentage of the overall density change from far upstream to far downstream (up to 
50 per cent in the example considered). 
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The Relation between Zero Errors and 
Bandwidths for Different Types of 
Rate-of-Turn Meters 


H. SUTCLIFFE 


(University of Bristol) 


Summary: In a conventional rate-of-turn measuring gyroscope the necessity 
to use a displacement transducer results in an appreciable zero error when 
a rapid response is required. Other types of instrument, employing alternating 
forces and force transducers, are also shown to have a relation between error 
and response time. An assessment is made of these relations and conclusions 
are reached concerning the relative merits of the different types. The results 
are given in terms of the estimated rate errors as functions of bandwidth. 


1. Introduction 


A rate-of-turn meter is an instrument which is mounted in a vehicle for the 
purpose of measuring the angular velocity of the vehicle about some axis. The 
instrument produces a signal which is normally used in connection with the 
mechanism which controls the vehicle. This present discussion is restricted to 
cases where the signal is proportional to the rate of turn. The signal is usually 
required to be electrical, and it follows that a rate-of-turn meter in its entirety must 
consist of all the sections shown in the schematic diagram of Fig. 1. The rate of 
turn is 2 and the output quantity of the instrument is v.,:. For an ideal instrument 
Vout=kQ, where k is a constant for the particular design. In a real instrument there 
will be four major deviations from this ideal state of affairs. Two of these, the 
variations in sensitivity due to variations in k and the phenomenon of saturation or 
non-linearity for large values of ©, will not be considered here. The two limitations 
to be discussed are the zero error of the instrument and the sluggishness of the 
response. 


The zero error is defined in terms of the indicated response to a constant 
applied value of 2. If vou=kA(2+.), then the rate error is ‘2. and is the 
indicated rate when °2 is zero. 


The sluggishness of the instrument is difficult to define in terms which are of 
general application to a variety of different instruments, and precision in the 


Received April 1957. 


May 1958 131 


H. SUTCLIFFE 


CORIOLIS 
MECHANISM 


FORCE PROPORTIONAL 


TRANSDUCER EITHER 
DIRECT OR VIA STRAIN 


wen ELECTRICAL SIGNAL 


ELECTRICAL NETWORK 


FicurE 1. The components of a rate-of-turn meter. 


definition must be sacrificed to make comparison possible. Assume that a test is 
made in which the input rate is given by 


cos qt. 


The output indication, neglecting any zero error, will be 


Vout = COs (qt 


When q is small, V..ax will be KQuax and @ will be zero. As q is increased it will 
pass through a value B at which the reduction of Vinax below the value KOQnax, or the 
value of ¢, or both, will be significant. This value of B will be taken to be the 
bandwidth of the instrument and the response of the instrument to a step function 
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The conventional rate gyroscope described in Section 2.1 is in common use 
and its limitations are well known to arise because the force produced is constant 
for a constant applied rate of turn. The transducer which converts this force into 
an electrical signal must therefore be of the displacement type used in connection 
with a mechanical strain’. Zero errors are associated with the transducer and 
bandwidth is restricted by the necessity to displace mass. These difficulties are 
removed when the force produced by the rate of turn is alternating, since piezo- 
electric transducers can be used which operate directly from the applied force. 
The principle has béen adopted in the design of several instruments')“’) and, if no 
other source of error had been found to occur in these, the conventional rate 
gyroscope would now be of historical interest only. 


Instruments which use an alternating force may be classified into two types, 
the dumb-bell gyroscope and the vibrating gyroscope. It is found in practice that 
in both these types, as in the conventional rate gyroscope, there is a rate error which 
is related to bandwidth. The nature of this relation is discussed in the following 
sections and an attempt is made to assess the relative merits of the three types on 


1 (t) will be similar to that of a single time-lagging network. That is to say the 
output will approximate to the shape vour=AQ2(1-e-"'). The response time may 
be defined as 1/B. The particular value of B to be allocated to different types of 
instrument is discussed later. 


a quantitive basis. 


All three types of instrument also have errors which are independent of band- 
width. A conventional rate gyroscope, for instance, will have errors due to mass 
unbalance and friction, while the other two types produce spurious forces at the 
transducer which are functions of the position of the moving parts. Errors arising 
from effects such as these become dominant only when a small bandwidth is 
adequate. A discussion of the situation under this condition is excluded from this 
paper, not because the subject lacks interest and importance, but because sufficient 
information about the relevant effects is not available. 


NOTATION 
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{radius of locus of point mass 
\amplitude of vibration of point mass 


distance from centre line 

bandwidth in radians /second 
moments of inertia about OX, etc. 
rotor moment of inertia 

sensitivity constant 

spring angular stiffness 


mass 
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B 
[,, etc. 
J 
k 
K 
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angular frequency of test signal 
torque about OY, etc. 
torque produced by error in instrument 
torque component which produces a signal 
Vou Output signal from instrument 
Vinax peak amplitude of electrical signal 
X,y,Z co-ordinates 
2 acceleration 
6, 6, etc. angle, angular misalignment. etc. 
phase angle of overall response 
angular frequency of spin or vibration 
rate of turn of instrument about OX, etc. 
error in measurement of 22 
The suffix ‘“‘max”’ denotes a peak value. 


Analysis of Errors and Bandwidths for Three Types of Instrument 
2.1. THE CONVENTIONAL RATE-MEASURING GYROSCOPE 


Figure 2 shows the instrument in one of the many possible mechanical arrange- 
ments. The rotor R has moment of inertia J and spin » about the OX-axis. The 
outer gimbal OG is fixed to the vehicle whose rate of turn 2, about the OY-axis 
is to be measured. Gyroscopic forces produce a torque T. about the OZ-axis, 
given by 

T.=JoQ,. 


The spring system § is displaced such that the angular displacement of the inner 
gimbal about OZ is 6, where @ is a small angle. Then if the angular stiffness of 
the spring system is K, the following relation holds. 


(1) 


The mechanism for transmitting the angle 4 will include a displacement 
transducer such as an E and J/-bar electromagnetic pick-off. A quantity 4, is 
postulated which is a measure of the misalignment of the instrument. It may be 
caused by imperfections in either the mechanical construction and assembly or in 
the transducer and its associated electrical circuits. The corresponding rate errof 
(), is given by 

/ (Jo). ‘ ‘ (2) 


The bandwidth of the response is limited by the moment of inertia 7 of the 
inner gimbal and rotor about OZ. The amplitude response is dependent on the 
amount of viscous damping, but a representative figure for the bandwidth can be 
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FiGurE 2. Conventional rate gyroscope. 


taken to be the frequency at which the phase of the response lags by 90°. This 
will be the case at the natural resonant frequency which is determined by spring 
stiffness and inertia only, and the bandwidth is given by 


It will be appreciated that a large value of K is required in order to provide a 

sufficiently large bandwidth, that is a stiffly sprung system must be used. It follows 

that for a given value of misalignment the rate error is correspondingly great. 

Quantitatively, by combining (2) and (3), 

The smallest possible ratio 7// would be attained with a disc rotor and a 

massless inner gimbal, when the value would be 1/2. A representative value which 
is obtainable in practice is unity, in which case equation (4) becomes 


2.2. THE DumB-BELL GYROSCOPE 


One possible form of this instrument is shown in Fig. 3. The rotor R spins 
with angular frequency » about the OZ-axis and may be considered, so far as the 
basic operation is concerned, to consist of two point masses m at radius a. The 
rotor is carried on an inner gimbal JG which is fixed in position by a pair of force 
transducers TT to an outer gimbal OG. This outer frame is fixed to the vehicle. 
The whole structure experiences rates of turn 2, and ©, about the OX- and OY-axes. 
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FiGureE 3. Dumb-bell gyroscope. Transducers TT measure torque about OY. 


Coriolis forces are transmitted between the vehicle and the spinning masses through 
the transducers, which are connected in such a way as to produce a signal which is 
a measure of the fluctuating component of the torque 7, about the OY-axis. 


The acceleration of the top mass in the OZ-direction is 2,, given by 


a, = 2aw [Q, cos wt + Q, sin wf]. 


Thus the total torque about the OY -axis is given by 


T,=2m2.4 COs wt 
=2ma’w (1 + cos 2f) + QO, sin 


The constant component is not measured by the transducers and the signal torque 
is given by 


The corresponding electrical signal from the transducers is amplified and 
passed to two demodulators (phase-sensitive rectifiers) in order to obtain two output 
signals proportional to ©, and 9,. The carriers which are also required at the 
demodulators must have the form cos 2wt and sin 2wt and can be derived from a 
pick-off operated by the angular position of the rotor. 


This analysis is concerned with magnitudes of error‘and it is convenient to note 
that , and ©, are components of an angular rate 2 about some line in the XY -plane. 
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CONTACT REGION 


FiGurE 4. One source of error in a dumb-bell gyroscope. 


Then, if Tax is the peak amplitude of the signal torque 7, in (6), the following simple 
result is obtained 


The instrument will have an error when a spurious signal coherent with the 
carrier frequency 2 is present in the signal amplifier. Such a signal may be 
caused by electrical interference and its elimination may present practical dijfficulties, 
but a more fundamental source of error is in the mechanical portion of the 
instrument. Imperfections in the rotor bearings can cause the rotor axis to deviate 
from the OZ-axis though small angles @, and 6, in the OX- and OY-planes. If 
6, or 6, has a component at 2w, then the torque 7, also has a component at 2. 
Figure 4 shows one of the more obvious sources of error. One of the rotor bearings 
is slightly elliptical, so that one end of the rotor axis shifts to and fro at twice the 
rotational frequency. 


Consider first the deviation 6,. Then, if the peak value of the component at 
2 iS max, the peak value of the angular acceleration of the rotor about OY is 
470, max. Let IJ, be a representative figure for the moment of inertia of the rotor 
about OY. Then the peak amplitude of the component at 2 of the torque is 
given by 
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Now consider the effect of a similar deviation 9,,.x. Torque about the 
OY-axis is now produced by gyroscopic action of the rotor, and the peak value at 
2 is given by 


where J is the moment of inertia of the rotor about OZ. 


Equations (8) and (9) ignore the fact that the rotor has a dumb-bell distribution 
of mass and are therefore simplifications of a somewhat more complex result. It 
will be realised, however, that in a practical instrument there will be a considerable 
amount of mass, with radial symmetry about the rotor axis, for which (8) and (9) 
are valid in an elementary manner. An example of a more detailed analysis is 
provided in the Appendix. The main purpose of this paper will be served by noting 
that, if the rotor structure approximates to a disc, then J in (9) will be nearly equal 
to 2/, in (8), and that a good approximation for /, is to take the value 2ma?. The 
phase of the error signal is of no significance, since the signals representing °, 
and {?, exist in the same circuit as vectors which are mutually at 90°. There is 
thus no possibility of eliminating errors by phase discrimination and it is reasonable 
to drop the suffixes and replace (8) and (9) by 


where 7. is the peak magnitude of the error signal. 
The corresponding rate error °2. is obtained by combining (10) and (7): — 


The available bandwidth is determined in practice by the presence of a signal 
in the transducer at frequency w, caused by rotor unbalance, by the desirability of 
arithmetical symmetry in the response of the signal amplifier about the frequency 
2w, and by the necessity to filter out the carrier component after demodulation. It 
has been found that B= /2 is representative of the highest practical value for the 
bandwidth. Then the equation which relates rate error, bandwidth and mechanical 
imperfections is 


2.3. THE VIBRATING GYROSCOPE 


Figure 5 shows the mechanical construction in one type of instrument, in 
which a pair of tines vibrates in the same way as a tuning fork. The situation is 
similar to the case of two masses oscillating sinusoidally along lines which in the 
ideal case pass through the axis OY about which torque is measured. Then if the 
amplitude of oscillation a is appreciably less than the distance b, the signal torque 
is given bv 


T, =2b6 x 2maw cos ; (13) 
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FicureE 5. One form of a vibrating gyroscope. 


This signal torque is converted by the transducer, amplified and demodulated in a 
phase-sensitive rectifier for which the carrier is coswt and is derived from tine 
velocity. The output of the demodulator is then proportional to ,,. 


One way in which an error torque can arise is shown in Fig. 5, where the 
lines of travel of the masses are displaced through angles 6, and 6, from their 
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correct orientation. The large acceleration of the masses along their lines of travel 
will produce an error torque 7., given by 


Tey=bmaw? (6,+6,)sinot, . . .  . (14) 


T., and T, differ in that they are sine and cosine functions of wt respectively. If, 
then, the control of phase angles throughout the whole instrument were perfect, the 
output signal after demodulation would contain no error. Phasing errors can arise, 
however, in every section of the instrument and their effect can be taken into 
account by adding a small angle @,, to the wt in (14). Then the component of 7., 
which is in phase with T, will have a magnitude 7. max given by 


T. max bmaw? (6, + 6.) . . . . (15) 


The rate error {2,, obtained from (13) and (15), is thus given by 


A more general summary, not restricted to any particular mechanical configura- 
tion, can be given, which leads to a similar conclusion. Any oscillating mass will 
have two accelerations, a wanted one due to the applied rate of turn, and an 
unwanted one due to its oscillation. The ratio of the amplitudes of these two 
accelerations, wanted to unwanted, is 20,/w. The two accelerations are each 
sinusoidal with time at angular frequency w, but in both relative direction and in 
relative phase angle they differ by 90°. The reduction of rate error depends on 
the success with which these two right angles can be maintained and, if they are in 
error by amounts @,,,, and ,,, the rate error is 

The suffixes ‘‘m’’ and “‘p” are reminders of the mechanical and phasing sources 
of error. 


The bandwidth of the instrument is restricted by the carrier frequency and the 
necessity for filtering after demodulation. A value B=w/2 will be taken as the 
largest bandwidth which will be practicable. Then the relation between rate error, 
bandwidth and imperfections is 


3. The Nature and Magnitude of the Sources of Error 


The three significant equations, (5), (12) and (18), are reproduced here for 
convenience. 


Equation of performance 


Conventional rate gyroscope 2, =6,B? 
Dumb-bell gyroscope 2,=96,8B 
Vibrating gyroscope 2,=9 


em ep 
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Each factor represents an error in alignment and it is necessary to define 
the nature of these errors before any conclusions can be drawn concerning the 
relative excellence of the three types, or of different designs of the same type. 
Suppose a batch of one design is being manufactured, tested and put into service. 
At the final stage of assembly, and before testing, each model will be aligned and 
adjusted to make {2, equal to zero, or very nearly so. The significant appearance 
of 6, and ©. will arise later when the model is subjected to changes of temperature, 
variations in “g,” wear of bearings, fatigue of highly strained metal, ageing of 
components and many other ill effects. A knowledge of @. and ©, for the design 
is thus derived partly from tests after assembly and partly from service trials. The 
errors have a correlation with the environment to which the instrument is subjected 
and the figure quoted for the design is that which can be achieved in the most 
trying environment. Any comparison between designs must be done in relation 
to the same environment. 


Within any one type it is clear that a smaller value of 6. and consequently of 
©, may be obtained for a particular design by a combination of ingenuity in 
conception, by the use of materials with suitable properties and by precision of 
workmanship. These requirements may be summed up in the word “cost.” It will 
be assumed here that different designs within one type will have the same value of 
6. if they have the same cost. 


To compare the relative excellence of different types, it is desirable to compare 


the performance of designs which have the same cost. The difficulties in the 
assessment of cost are obvious and it must be emphasised that the quantitative 
comparison which follows is based entirely on the author’s restricted experience 
and intuition. The figures for rate errors may well be incorrect by a large factor 
for the dumb-bell and vibrating gyroscope, since the techniques involved in their 
design are in their infancy and conclusions must be drawn from one or two samples. 
This proviso, however, does not detract from the significance of the way in which 
rate error varies with bandwidth in the different types. The quantitative comparison 
relies on the use of arguments such as the following. 


The thermal coefficient of expansion of most metals''”’ is in the region of 10~° 
per °C. Thus a structure of uniform material with a square cross section, with 
one pair of opposite faces differing in temperature by 5°C, will experience a tilting 
of the relative angle of the other two faces through an angle 1/20,000. In a similar 
shaped structure whose opposite sides are made of materials whose coefficients of 
expansion differ by 10-* per °C, the same tilt is reached for a temperature change, 
uniform in the whole structure, of 50°C. 


The construction of two separate parts of a mechanism to limits of half of one 
ten thousandth of an inch, at a radius of one inch, represents very precise work- 
manship. The corresponding angle is 1/20,000. 


It may be expected then that the maintenance of the alignment of angles to 
better than 1/20,000 in any instrument will be achieved only by a considerable 
increase in cost. This will be particularly true for the vibrating gyroscope, where 
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the error angle @.,, arises in a part of the mechanism which must be free to vibrate. 
At equivalent cost, then, a prediction for 6, for the conventional gyroscope is 
1/20,000 and the value 4.,, for the vibrating gyroscope may be expected to be 
somewhat larger. For the dumb-bell gyroscope a prediction of 6, is rather difficult, 
since there is no general information about the second harmonic motion of high- 
speed bearings. A very approximate estimation can be made by assuming that 
the running clearance in the bearings will be in the region of 0-0001 in. Then, 
with bearings two inches apart, the rotor axis cannot change through an angle 


greater than 1/10,000. Thus 4 cannot be greater than 2/= times this and is likely 
to be appreciably smaller. 


Turning now to specific instruments, from which conclusions can be drawn, 
a conventional rate gyroscope is considered first. A typical design has a wheel 
speed of 12,000 r.p.m. and a resonant undamped frequency of 100 c./s. It passes 
a severe environmental test for which a rate error of 1° per second is allowed. 
The corresponding value of 6, is 1 / 18,000. 


Model dumb-bell gyroscopes have been built and tested. Values of 6 were 
measured both directly and by inference from rate errors. Good correlation was 
obtained over a range of wheel speeds and a figure of 3:5 x 10-* was arrived at for 
6. for the same conditions of cost and environment as in the previous case. 


Experimental vibrating gyroscopes have been built and tested. A typical value 
of rate error was found to be 2° per second when the vibrational frequency was 
100° c./s., so the value of (6.,4.,) in the models was in the region of 1/9,000. The 
models, however, had a low cost in the sense which has been defined, the sources of 
error were well established, and it is believed that designs whose cost was equivalent 
to the other two examples would have values of (6.,4.,) in the region of 1/30,000. 


The result of these considerations is to suggest the following values for the error 
quantities in equations (5), (12) and (18). 


Conventional rate gyroscope: 6.=1/20,000 
Dumb-bell gyroscope: 1 /300,000 
Vibrating gyroscope: FomFep = 1 / 30,000. 


It is now possible to make a direct comparison between the three types of 
instrument. Figure 6 shows the minimum rate error available in the three types, 
assuming equivalent cost and environment, plotted against bandwidth. Logarithmic 
scales are convenient for the presentation. For the conventional rate gyroscope a 
rotor speed of 24,000 r.p.m. (400 c./s.) is assumed whatever bandwidth is required, 
whereas in the other two cases it is assumed that rotor speed or vibrational frequency 
is as low as is consistent with obtaining the required bandwidth. The highest rotor 
speed for the dumb-bell gyroscope is taken to be 24,000 r.p.m., which implies an 
upper limit of 200 c./s. for the bandwidth. 
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BANDWIDTH IN C/S. 


FicurE 6. Relation between rate errors and bandwidth in comparable designs of the three 
types of instrument. 


4. Conclusions 


The first point arising from the results shown in Fig. 6 is that the dumb-bell 
gyroscope appears to be outclassed by the conventional type, although the situa- 
tion could be modified by a fortunate discovery in the subject of elimination of 
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second harmonic play in high-speed bearings. The dumb-bell gyroscope would 
then be the best type for applications which require bandwidths up to the region 
of half the upper limit of wheel frequency. 


At present the two main contestants are the conventional and vibrating types 
of instrument. The conventional type is superior for bandwidths up to the region of 
300 c./s., while the vibrating type would be superior if wider bandwidths were 
required The increase of error with bandwidth is an inexorable feature of all 
the instruments. 


The results become invalid for small bandwidths. Fig. 6 suggests that the zero 
error in all three types decreases continuously, whereas each type of instrument 
has a zero error, independent of bandwidth, which becomes dominant at small 
bandwidths. With regard to the conventional rate gyroscope, quantitative informa- 
tion about the performance under these conditions is in existence. An investigation 
of the performance of the other two types, using models in which the moving parts 
had low rotational or vibrating frequencies, could yield some interesting and 
perhaps valuable results. 
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Appendix 
ANGULAR VIBRATION OF THE DUMB-BELL GYROSCOPE AT TWICE ROTOR SPEED 


The practical case will be considered in which there is a concentration of mass 
towards the dumb-bell axis, instead of the simpler case of two point masses. Suppose 
that mass is concentrated towards an axis OX, referred to the rotor. The rotor will be 
assumed to spin about an axis OZ, and to have symmetry about the X,Z,- and 
Y,Z,-planes. Let the rotor spin at'angular frequency w with respect to axes OX, OY. 
OZ referred to the frame or gimbal, such that OZ coincides with OZ, at all times and 
that x=x,, y=y, at t=0. Then the position and velocity of an element of mass 6:7, 
with respect to the frame, is given by 


£=2, 
X=X, COSwf-y,sinwl, x =o[-x, sin wt - y, cos wf] 


y=X, SiN wf+y,coswt, y =w [x, cos wt y, sin wf]. 


The transducer is constructed so that it measures the torque 7, about the OY-axis. 


(i) The Effects of Rates of Turn Q, and Q, 


T, = (20, x) Ldmx (20, y). 
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Now 3émx,y, is zero, by the symmetry of construction of the rotor. 


Therefore T,= (x,? - sin Zot (x,? + - 
(x,? — cos 2ut. (19) 
Let J be the rotor moment of inertia about OZ, and let /,, the effective dumb-bell, be 
dm 


Then T,= - wl, [O, cos 2wt +Q, sin 2wt]-JoQ,. (20) 


Consider an oscillation of the rotor axis through an angle 6, at an angular frequency 
2m such that 0,=200, max Sin 2wt. Then the component of 7, at 2 will have an 
amplitude given by 


(ii) The Effect of Angular Acceleration about OY 
Let the angular acceleration be 0, 
Then T, (2? 
=1,38m [z,? + cos? wt + y,? sin? wf] 


= [4(z,? + X02) + + Yo?)] + 40, 38m (x,? — C08 (22) 


Then, if the moments of inertia of the rotor about its axes OX, and OY, are /,,, and 
amplitude given by 
I,. respectively, 


Consider an oscillation about OY of peak value @, ,,,, such that 0, = 406, max SIN 2wf. 
The component of 7, at 2w will have an amplitude given by 


Tes max — 2 max’ be (24) 


(iii) Rate Errors 


Let rate errors (),, and {2,, correspond to torques T,, and T,,. 


Then = 2086 and = max no + Lae, 


__ In the ideal case there will be mass only along OX,, and both J/7, and (1, +1,.)/Ta 
will be unity. They cannot be smaller. Then the limiting value of rate error is 


209. max» A practical case had values as follows :— 
§=2325, 1680, 7/0 gmem-, 
from which 1,=910 gm. cm.? 
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Then J /1q=2-55, + = 2°69. 


It will be seen that the numerical factor 4 in equation (11) is representative of what 
can be achieved in practice. 


It will be seen from (23) that any angular acceleration of the vehicle appears as a 


rate error. This feature is common to all three types of instrument. 
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Photoelastic Investigation on Plates with Single 
Interference-Fit Pins with Load Applied 
(a) to Pin only and (b) to Pin and Plate 

Simultaneously 


H. T. JESSOP, C. SNELL and G. S. HOLISTER 


(Department of Civil and Municipal Engineering, University College, London) 


SumMMARY: The reduction in stress concentration factors resulting from the 
introduction of an interference-fit pin in a circular hole in a flat plate to which 
a simple tension is applied has already been recorded. The same result is 
found to occur to an equally marked extent, and qualitatively in a very similar 
way, when the load is applied to the pin itself. When varying loads are applied 
to the pin or the plate, the effect of the interference is to produce a rise in the 
mean stress level at critical points on the hole boundary with a marked fall 
in the oscillatory stress. There appears to be an optimum value for the 
“interference stress” for given applied loading conditions at which the endurance 
of the plate will have its highest value. 


When loads are applied simultaneously to pin and plate the critical stress 
on the hole boundary is found to be the mean of the stresses which would be 
produced by the total load applied separately to the pin and to the plate. 


1. Introduction 


This is the fourth stage of the photoelastic tests in the programme of research 
initiated by the Royal Aeronautical Society and outlined in Ref. 1, and it concludes 
the investigation of the “single hole” case. Previous results have been published in 
Refs, 2, 3 and 4. This investigation deals with plates with circular holes filled by 
interference-fit pins with loads applied (a) to the pin only, and (b) to pin and plate 
simultaneously. 


In this investigation a series of tests was first carried out on plates with loads 
applied to the pin only, the intention being at that time to make a few “spot” 
observations with loads on pin and plate simultaneously. It was decided, however, 
to carry out in the second stage a series of tests covering the whole range of com- 
bined loads. In this series, the previous observations under pin load only and 
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plate load only were repeated, and were found to give very good agreement with 
those of the earlier tests. In this report the results are all taken from the combined 
load observations. 


NOTATION 
E Young’s modulus of plate 
E, Young’s modulus of pin 
D_ width of plate 
d,d’ diameters of hole and pin respectively 
(d’-d)/d_ degree of interference 
t mean tensile stress in undrilled plate due to applied tension 


t’ nominal mean tensile stress across section through hole centre, due to 
applied tension, = Dt/(D- d) 


tm = maximum tensile stress on hole boundary 


Sm= Maximum shear stress on hole boundary 


The stress concentration factors derived here are given in terms of the increase 
in peak tensile or shear stress produced by the applied tension. 


Thus : — 
S.C.F. for tensile stresses =5t,,/t’ 


S.C.F. for shear stresses =3s,,/f 


2. Conditions of Tests 
As in previous tests, the models were Araldite bars, 1 in. wide and } in. thick. 


The pin load was applied through two thin steel strips which passed over the 
projecting ends of the pin and were attached to a skeleton box frame. These strips 
can be seen in the photographs of Fig. 1. The ends of the pin were turned down, 
in order that the strips should obscure as little of the hole boundary as possible. 


The tests in this series were confined to plates with two sizes of hole, namely 
} in. and ? in. diameter—these representing the ratios of hole diameter to plate 
width most frequently used in practice. The pins used throughout the test were 
of Bakelite B.166, a paper base Bakelite giving a modular ratio E,/E of 2-15, the 
nearest obtainable to the 3/1 ratio required to represent a steel pin in a light alloy 
plate. (In previous tests no significant difference in stresses was observed for 
different ratios of E,/E). 


Three degrees of fit were investigated; push fit and two degrees of interference. 
namely 0-003 and 0-006, corresponding approximately to maximum elastic 
shear stresses in a light alloy bar of 20,000 lb./in.* and 40.000 Ib./in.? respectively. 
These are specified as “low” and “high” interference respectively. 
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(c) (d) 


FiGurE 1. Typical fringe patterns in a plate with a ? in. diameter pin. The actual hole boundary 
is masked by the pin. The fringe numbers indicated are those on the boundary obtained by 
extrapolation along a normal to the boundary, the scribed 4 in. diameter circle giving reference 
points for measurements. 

(a) Push-fit, load 16 x 10* Ib. /in.? 

(b) Push-fit, load 32 x 10° Ib./in.? 

(c) Low interference (0:003), load 16 x 10* Ib./in.? 

(d) High interference (0-006), load 32 x 10* Ib./in.? 
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For each degree of interference, a series of progressively increasing pin loads 
was applied, firstly with zero plate load and then with superimposed plate loads of 
three different magnitudes. 


For each loading and pin-fit condition, the stresses at three points on the hole 
boundary were analysed. These points were: (a) on the horizontal diameter of the 
hole; (b) on the vertical diameter directly below the pin load, and (c) at the lower 
end of a diameter inclined at approximately 45° to the vertical. 


The results of these tests are given in Tables to V. The loads applied to the 
photoelastic model have been converted in these tables to loads which would 
produce the same strains in a light alloy plate, and the stresses in the model have 
been scaled up in proportion. Thus the tables give loads and stresses for a light 
alloy plate in which the degree of interference in the pin is the same as that in the 
photoelastic model. 


3. The Stresses due to Pin Load Only 
3.1. THE STRESS DISTRIBUTION AROUND THE HOLE BOUNDARY 


Figure 1 shows photographs of typical fringe patterns, while Figs. 2 and 3 
show respectively the distributions of shear stress and tensile stress around the hole 
boundary for two different loads at each degree of interference. The only 
separations of stresses carried out for the upper half of the hole boundary were 
those for the point on the vertical diameter in the cases of push-fit and high 
interference. The tensile stresses indicated by the dotted broken curves around 
the upper half of the hole boundary in Fig. 3, therefore, are only approximate, 
except for the four points indicated in Figs. 3(a) and 3(c). In the push-fit case lift 
from the pin occurs and the boundary stress (a compression one) is given directly 
by the fringe pattern. An examination of these figures gives the following results. 


(i) In the push-fit case the maximum shear stress occurs at a point about 30° 
below the horizontal diameter, while the tensile stress has a much more pronounced 
maximum on the horizontal diameter. Thus the conditions at the point of maximum 
shear stress under the smaller of the two loads illustrated are 


Shear stress 25 x 10° Ib./in.? 

Circumferential tensile stress 37 x 10° Ib./in.* 

Radial compressive stress 13 x 10° Ib./in.? 
while at the horizontal diameter they are 

Shear stress 23 x 10° 1b. /in.? 

Circumferential tensile stress 46 x 10° Ib./in.? 

Radial stress (approximately) 0 Ib. /in.? 


Since in practice such a loaded lug never fails at a point below the horizontal 
section, it would appear that the small (8 per cent) decrease in shear stress at that 
section is more than offset by the 25 per cent increase in tensile stress. 
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(ii) With low interference and high load the stress distribution characteristics 
are similar to those in the push-fit case, but under low load both tensile and shear 
stresses tend to maintain almost constant values around a large section of the hole 
boundary. the maximum shear stresses having nearly the same values as in the 
push-fit case (see Figs. 1(a) and 1(c)), while the tensile stresses are considerably 
lower. With the two loads illustrated, the circumferential stresses remain tensile 
all round the boundary. 


(iii) With high interference under both high and low loads the stresses, both 
shear and tensile, are very nearly constant around the greater part of the circum- 
ference, only falling towards the top of the hole. At the extremity of the horizontal 
diameter the shear stresses are somewhat higher than in the push-fit case, but the 
tensile stresses are lower, the differences in each case being greater under the lower 
load. The relative magnitudes of the shear stresses are shown by the fringe numbers 
in Figs. 1(b) and 1(d). 


3.2. THE EFFECT OF INTERFERENCE UPON THE SHEAR STRESSES 


Tables I, II and III give the shear stresses at the three points investigated for 
the whole range of loading for the two pin sizes. The results for pin load only have 
been taken from the appropriate columns of these tables. 


It may be noted that, in this case, the circumferential stress becomes 
compressive at both the top and bottom of the hole. At the top of the hole “lift” 
of the plate from the pin will occur at all loads and there will be regions of very 
small radial contact pressure in the upper half of the periphery at which slip and 
consequent fretting will occur. 


TABLE I 
SHEAR STRESSES ON HORIZONTAL SECTION 
(Unit=10° 1b./in.*) 
} in. Diameter Hole 


Push-Fit ‘Pin i Low Interference (0-003) High Interference (0-006) 
P at Plate Load Plate Load Plate Load 
Loa 0 133 267 400 | O 133 26:7 40:0 0 133 26:7 40:0 
0 0 17 34 SO 2% 36 40 41 63 
G7 | 13 2 44 62 21 33 | 66 40 44 56 vi 
133 | 26 40 55 75 28 «45 65 80 50 65 90 
20:0 40 50 69 90 39s «60 78 96 44 58 78 100 


267 | 53 60 85 105 54. 78 90 110 50. 70 95 110 


3 in. Diameter Hole 


Pin Plate Load Plate Load Plate Load 
Load | 9% 160 320 480) 0 160 32:0 480)! 0 160 32:0 48-0 
0 | 0 19 38 55 20 25 40 50 40 42 50 57 


8-0 | its 6°30 50 66 21 31 47 65 41 a4 55 67 
16:0 | 23 43 60 78 25 41 62 80 | 43 46 60 77 
24:0 | 34 55 72 89 =| (34 55 75 94 | 47 54 70 88 
32:0 | 46 65 84 100 48 68 88 106 | 53 68 80 100 
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cs TABLE II 
ar SHEAR STRESSES ON 45° SECTION 
(Unit=10° Ib./in.2) 
he 1 in. Diameter Hole 
ly Push-Fit Pin Low Interference (0-003) | High Interference (0-006) 
ile Pin 4 Plate Load | Plate Load Plate Load 
Load | 9 133 267 400 0 133 267 400 | 0 133 267 400 
th 0 | 0 7 4 2 2020 «#242 «2 | 4 39 «#3440 «(41 
n- 67 | 15 19 23 30 2s 622 35 40 | 42 40 43 48 
| Iss | 30 SI 32 40 ; 30 31 48 53 | 45 45 48 55 
la 20:0 45 43 43 48 | 38 47 61 65 50 52 55 64 
he 267 | 60 50 55 55 49 60 70 £78 52 62 62° 73 
TS 3 in, Diameter Hole 
Load | 9 160 320 480) 0 160 320 480) 0 160 32:0 480 
0 0 6 12 18 20 20 20 28 40 40 40 40 
or 80 | 11 18 20 26 | 22 22 22 34 | 42 41 41 43 
ve 16:0 | 23 25 28 35 | 26 25 30 40 45 44 43 48 
24:0 35 35 38 43 33 34 40 45 47 46 46 54 
| 32:0 | 46 47 48 52 43 45 50 51 50 50 50 60 
ies 
ry 
nd 
TABLE III 
SHEAR STRESSES ON VERTICAL SECTION 
(Unit=10° 1b. /in.?) 
4 in. Diameter Hole 
Push-Fit Pin Low Interference (0-003) | High Interference (0-006) 
) rol | Plate Load Plate Load Plate Load 
7 0 133 267 400 0 133 267 400 0 133 267 400 
) 0 | 0 -6 -13 -19 20 15 8 4 40 34 25 18 
o7 | 3 2 1 22 17 19 14 43 35 30 23 
13-3 18 12 13 12 27 24 30 24 46 40 36 30 
22 18 33 36 39 33 50 47 43 39 
26°7 | 36 28 30 24 41 55 48 42 Sz 6§5 52 51 
in. Diameter Hole 
Pin Plate Load Plate Load | Plate Load 
0 Load | 9 160 320 480 0 160 320 480) 0O 160 320 48-0 
0 |; 0 -5 -10 —-—15 20 12 7 4 40 35 25 16 
80 | 6 4 2 0 21 15 9 8 42 36 y4 | 18 
16:0 | 11 12 9 8 22 19 14 13 | 45 37 30 22 


240 | 17 20 14 14 | 25 24 21 17 | 48 40 33 26 
320 | 23 26 20 20 30 29 26 22 =| «+50 43 37 30 
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(Q) HORIZONTAL SECTION 
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Ficure 4. Shear stress concentration factors under pin load, (a) on the horizontal section. 
(b) on the 45° section, for 4 in. and ? in. diameter holes. 


As for the plate load“, the effect of the interference is found to be directly 
related to the ratio of interference stress to applied tension, and can best be 
illustrated by graphs of the stress concentration factor plotted against this ratio. 
Figure 4 shows such graphs for the shear stresses, (a) on the horizontal section, and 
(b) on the 45° section. 


The reduction in S.C.F. is, of course, offset to some extent by the initial 
interference stresses, but the net result may be shown on a graph of maximum total 
stress as a function of applied tension. Figure 5 shows this graph for the } in. 
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FicurE 5. Maximum and minimum shear stresses on the horizontal section when the applied 
pin load alternates between zero and a given value #’. Hole diameter 4 in. 
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FicurE 6. Peak shear stresses on the horizontal section for various degrees of interference 

under pin load. Hole diameter 4 in, The broken line denotes the estimated limits of load for 

which the stress concentration factors under reversed load will have the same values as for 
positive load. (See Section 5 (iii).) 

(The interference figures given are those for a steel pin in a light alloy plate. For a light alloy 

pin the interference to give the same initial stress will be approximately 4/3 of that in a steel pin) 
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diameter hole. The straight line in this figure represents the initial interference 
stress, and the graph is designed to show the mean stress level and the amplitude of 
stress oscillations when the load on the pin alternates between zero and a given 
fixed value. 


For the general case in which there is a mean tensile load on which an 
alternating load is superimposed, a more useful way of presenting the results is 
illustrated in Fig. 6. Using one of the S.C.F. curves of Fig. 4, curves of peak stress 
against applied tension are deduced for various values of interference and, from 
these, it is a simple matter to read the mean value and the amplitude of the stress 
for any given loading conditions. 


TABLE IV 
TENSILE STRESSES ON HORIZONTAL SECTION 
(Unit = 10° Ib. /in.) 


1 in. Diameter Hole 


Push-Fit Pin | Low Interference (0-003) High Interference (0-006) 
pose Plate Load | Plate Load Plate Load 
13-3 267 | 13-3 26°7 0 13-3 26-7 
0 0 26 yi ee 32 49 | 44 46 51 
53 75 32 70 39 60 
26°7 106 — 145 81 _— 160 36 _ 150 
2 in. Diameter Hole 
Pin | Plate Load | Plate Load | Plate Load 
Load| © 160 320 | 0 16:0 320 | 0 16-0 32:0 
0 | oO 26 52 | 23 27 40 | 48 48 51 
160 | 46 74 — | 27 57 — | 4 54 — 
92 150 67 — 115 a7 120 
FI 
TABLE V 
TENSILE STRESSES ON 45° SECTION 3 
(Unit = 10" Ib./in.?) 
1 in. Diameter Hole 
Push-Fit Pin Low Interference (0°003) | High Interference (0-006). 
Plate Load Plate Load | Plate Load 
13-3 26°7 | 0 13-3 26°7 0 13:3 26:7 
0 0 3 6 | 2 20 19 | 40 40 40 e 
i373 | 2 15 28 20 — | 4 40 
26'7 56 _ 30 39 — 35 41 —_ 40 1 
: in, Diameter Hole h 
Pin Plate Load | Plate Load Load 
0 160 320 | 0 16:0 32:0 0 16-0 32-0 
0 0 2 s | 20 20 20 40 40 ~— 40 A 
16:0 21 15 — | 29 20 _ 43 40 _— b 
32:0 42 — 30 36 — 30 53 -- 40 
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FiGURE 7. Tensile stress concentration factors under pin load, (a) on the horizontal section, 
(b) on the 45° section, for 4 in. and 3 in. diameter holes. 


3-0 


3.3. THE EFFECT OF INTERFERENCE ON THE TENSILE STRESSES 


Separation of the stresses was carried out for a limited number of load 
conditions and the tensile stresses obtained on the horizontal and 45° sections are 
given in Tables IV and V. 


The accuracy of these results is not as high as that of the shear stresses, 
especially for the 45° section, where the isoclinic fringes in the photoelastic model 
were frequently not well defined. The scatter in the S.C.F. graphs of Figs. 7(a) and 
7(b) indicates, however, maximum divergences of the order of 15 per cent of the 
highest stress. 


Figure 8 shows the graphs of maximum and minimum tensile stress on the 
horizontal section for a } in. diameter hole under application and removal of load. 
A reversal of sign of the S.C.F. for a high ratio of interference to load is indicated 
by the crossing over of the two curves. This reversal, which indicates a reduction 
in tensile stress on the boundary with the addition of a tensile load, might appear 
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FiGuRE 8. Maximum and minimum tensile stresses on the horizontal section when the applied 
pin load alternates between zero and a given value r. Hole diameter in. 


at first sight to be impossible, but if we postulate an interference great enough to 
prevent relative motion of the pin and plate surfaces, then it can be seen that either 
a small vertical compressive strain in the pin at this point or a small increase in the 
normal pressure between the surfaces will require a vertical compressive strain in 
the plate to maintain the contact conditions. 


4. The Stresses under Combined Pin and Plate Loads 


The figures of Tables I, II, [1V and V, deduced from the individual photoelastic 
observations, were first analysed to check their consistency. The stresses were 
plotted against load (a) for varying pin loads at each plate load, and (b) for varying 
plate loads at each pin load. Smooth curves were then drawn which gave the best 
“fits” to the plotted points, while keeping the readings of corresponding points on 
the two sets of curves identical. 


All the shear stresses on the horizontal section gave equally good curves, with 
divergencies of the same order as in Fig. 9. The divergencies of the shear stresses 
on the 45° section and the tensile stresses on both sections were somewhat greater 
than these. 


Figure 10 shows maximum shear stress values on the horizontal section for 
various proportions of pin and plate load. The greatest number of observations 
was available for the case of equal pin and plate loads, and these are found to lie 
sensibly mid-way between the pin and plate load curves. To the same degree of 
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FIGURE 9. Graphs of shear stress on the horizontal section under increasing pin and plate loads, 
showing the small adjustments needed to produce smooth curves giving consistent values of the 

stress under all load conditions. 

(Taken from Table I for 4 in. diameter hole and 0-006 interference) 
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FicurE 10. Values of maximum shear stress on the horizontal section for different proportions 
of pin and plate load. Hole diameter } in. 
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accuracy the other observations divide the range in the inverse ratio of the propor- 
tions of pin to plate load. This proportional sharing of the total stress is found 
for both shear and tensile stresses on the horizontal section. Thus, under equal pin 
and plate loads ?, the stresses with any degree of interference will be the means of 
those for pin and plate loads each of 2/ with the same interference. 


The stresses on the 45° section do not follow this pattern. The curve for equal 
pin and plate loads lies much nearer to the plate load curve, and this effect is even 
more marked in the tensile stresses. Thus the effect of superposing a plate load is 
to reduce the effective S.C.F. of the pin load for the 45° section. 


5. Applicability of Results to Practical Problems 


(i) All the results obtained relate to an elastic distribution of stress. In 
cases where the elastic limit of the material is exceeded, the peak stresses indicated 
in these tests may not be reached, but the assumption is made that endurance under 
fatigue is dependent upon the peak stresses which would have been reached if local 
plastic yield had not occurred. 


(ii) No account is taken of the effect of fretting, which may be an important 
factor, especially in cases of low interference-load ratios. 


(iii) The results are not necessarily directly applicable in cases of reversal of 
sign in the applied load. It seems probable that, so long as the interference-load 
ratio is high, the contact conditions will not alter seriously with reversed loading 
and there will always be a reserve of radial compression allowing of the transmission 
of tensions across the hole boundary. In that case the stress-concentration factors 
under reversed loading should be the same, leading to load stresses of the same 
magnitude but opposite sign. For low interference-load ratios, however, where 
“lift” from the pin tends to occur under tension, the S.C.F.’s under reversed loads 
will be entirely different. 


Thus in the graphs of peak stress shown in Fig. 6, the curves may be continued 
beyond the vertical axis in inverted form for a limited distance. The broken line 
represents an estimate of the range within which the interference-load ratio is high 
enough to justify this. 


(iv) Within the limitations imposed by the factors mentioned in (ii) and (iii), 
the results obtained should help the designer to estimate the endurance of joints 
containing interference-fit pins, or to determine the degree of interference which 
will give maximum endurance under given loading conditions. 


(v) It is hoped also that these results will be of value in investigations into 
the basic causes of fatigue failure at bolt and rivet holes. With much more precise 
information than has hitherto been available on the peak stresses in such features, 
it should be possible to devise fatigue tests to explore systematically the separate 
effects of the various factors affecting fatigue failure. 
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Appendix 


Another method of presenting the results on interference fits has been proposed 
by H. L. Cox. In this method, charts are drawn which show the ratios of mean stress 
and stress range with different degrees of interference to the corresponding stresses 
with a push fit. These charts enable the designer to predict the fatigue endurance of 
interference-fit joints from fatigue data on similar joints with push-fit pins. 


Separate charts are needed for plate load and pin load, and each chart applies 
strictly only to a given ratio of mean load to load range. The variation of the stress 
ratios with the load ratios, however, is very small over a considerable range, and charts 
for three different values of load ratio should be sufficient to cover the range of loading 
conditions normally occurring in aircraft practice. 


Two such charts are shown in Figs. 11 and 12. These were prepared at the 
National Physical Laboratory, Ph.(E) Division, by A. F. C. Brown and A. G. Cole, who 
also carried out the additional photoelastic tests to provide data for the + in. and 4 in. 
diameter holes under pin load. The authors are indebted to Messrs. Cox, Brown and 
Cole for permission to include the charts in this paper. 


The notation used in the charts is as follows : — 
t plate thickness (in.) 

D plate width (in.) 
d_ pin diameter (in.) 
P tensile load applied to plate (Ib.) 

P.nax Maximum tensile load applied to plate (1b.) 

P.xin Minimum tensile load applied to plate (Ib.) 
M_ mean.tensile load applied to plate, = 4 (Pmax + Pmin) (Ib.) 
R_ stress ratio, =(Pin/Pmax) 


f average tensile stress across net area of plate due to mean tensile load, 
=M/{(D —d) t} (Ib./in.”) 
q Shear stress at transverse boundary of hole caused by applied load P 
with push-fit pin (Ib. /in.*) 
q shear stress at transverse boundary of hole caused by applied load P 
with interference-fit pin (Ib. /in.*) 
q, shear stress at transverse boundary of hole caused by interference only 
(Ib. /in.?) 
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R=0°5 
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05 Fe) 15 20 25 30 35 40 


Ficure 11. Ratios of mean stress and stress range on the horizontal section with different 
degrees of interference to corresponding stresses with a push-fit pin. under an oscillatory plate 


load whose minimum value is 0°5 of its maximum. 
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Load applied to Plate only). The Aeronautical Quarterly, August 1955. 


4. Jessop, H. T., SNELL, C. and HouisTer, G. S. Photoelastic Investigations on Plates with 
Single Interference-Fit Pins with Load Applied to Plate Only. The Aeronautical Quarterly, 
November 1956. 
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The Shock Pattern of a Wing-Body -~ 


Combination, Far from the Flight Path |"*' 


F. WALKDEN 


In ec 
(The University, Manchester) 
and | 
dista 
SUMMARY: The position and strength of the front shock wave at large distu 
distances from a wing-body combination, are deduced from the linear duce 


theory for the combination, using a method developed by Whitham. The 
combination consists of a body of revolution and a wing which has as fc 
thickness and is lifting. The effects of interference between the flow over 
the body and the flow over the wing are included. In any direction the flow 
far from the wing-body combination is equivalent to the flow past a body of 


revolution determined from the configuration of the combination. The the s 
modified formulae for unsteady flow are given and some results are evaluated 
for the combination of a body of revolution and a delta wing with subsonic 
leading edges. proc 
treat 
1. Introduction equi 
Tevo 
This paper is a contribution to the theory of the shocks produced by aircraft all ¢ 
and missiles in supersonic flight. In recent papers by Whitham“:” a method has 
been developed for deducing this information in certain cases from the linear theory 
of the disturbance produced in the air. In Whitham’s first paper this method is is th 
applied in detail to a body of revolution in uniform motion, and the second paper of tl 
includes the extension to a thin symmetrical wing in uniform motion. Rao has dist 
shown how the method can be extended to deal with the non-uniform motion of a 
body of revolution. inve 
Here it is intended to investigate the flow far from a general wing-body com- tion: 
bination which is moving with supersonic speed. It is assumed that the body is 
axi-symmetric, with the main stream direction as axis of symmetry; the wing may je 
be inclined at a small angle to the main stream. In the first instance uniform motion x 
is assumed, so that the problem can be treated in the steady state. However, the - y 
further extension used by Rao to obtain the results for non-uniform motion from a 
those for steady motion can still be applied. For the sake of simplicity, attention 
is concentrated on the flow at large distances, but it should be noted that for some 
body shapes the theory gives results at all distances from the body surface. 4 invc 
Received March 1957. p bec 
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For the steady-flow problem. cylindrical polar co-ordinates (r, 9, x) are 
introduced, with origin at the intersection of the leading edge of the wing and the 
axis of symmetry of the body; the axis of symmetry of the body is also the 
axis of x. 


Let U.,¢ be the velocity potential of the disturbance. From the linear theory 
it is found that, far from the body, # takes the approximate form 


_ f(x—Br, 6) 
(1) 


In equation (1), B= /(M,,* - 1), where M,, is the Mach number of the main stream. 


This potential is made up of contributions from the body alone, the wing alone, 
and the interference potential. All these contributions take the same form at large 
distances. This is in accordance with the theory of geometrical acoustics; as the 
disturbance spreads out from the body the total energy flux remains constant and, 
since the area is increasing proportionally to r, the amplitude falls off like 1/ /r. 
as for cylindrical waves. 


The function f gives the initial profile of these waves and it is determined by 
the shape of the body. 


The method of finding the shocks from equation (1) follows the standard 
procedure explained by Whitham. Briefly, the flow in each plane 6=constant is 
treated separately; the flow in such a plane may be interpreted as being due to an 
equivalent body of revolution and the theory can be developed as for a body of 
revolution. A more involved treatment is necessary to determine the shocks at 
all distances; the simplification applies only at large distances. 


Rao™ verified that, for a body of revolution, the function f for unsteady motion 
is the same as the function for uniform motion. The effect of non-uniform motion 
of the body appears in the amplitude factor, which governs the way in which the 
disturbance decays with distance from the body. 


Whitham found the dependence of f on 6 for a thin symmetrical wing. The 


_ Investigation will be completed here by finding f for general wing-body combina- 
_ tions and then applying it to a study of the shock waves. 


For a body of revolution f is determined as a certain integral involving S (x), 
the area distribution of the body. For the symmetrical wing f depends in the same 
Way upon a function S (x), which is found by integrating the thickness distribution 
of the wing plan form; the choice of lines depends upon the direction being 
considered. 


__ For the anti-symmetrical problem, it will be shown that f is a similar integral 
involving the lift distribution over the wing plan form. The lift distribution can 
be obtained in terms of the wing shape by standard methods™. 
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For a winged body the perturbation potential in the linear theory is the sum 
of the potentials of the wing and body considered separately, together with an 
interference potential; this interference potential arises from the effects of the 
interaction between the two flows. 


Chester has studied the problem of interference. At large distances the 
main contribution is from a term which is equivalent to a distribution of sources 
along the axis of the body; the strength of this source distribution is proportional 
to the mean of 09,/dr taken over the body surface. U.y is the perturbation 
potential of a two-dimensional wing whose profile is equal to the profile of the 
centre section of the given wing. The derivation of this term is different from 
Chester’s method and is included in Section 2. 


In a private communication, Warren has considered the effect of lift on shock 
waves. He represents the effect of lift as that of a distribution of doublets along 
the axis; he chooses the strength of the doublet distribution to get a realistic total 
lift (although this leads to an impossibly large incicence for the body of revolution 
which is his model) and obtains plausible results for the shock wave behaviour. 


The method used in this paper gives the same kind of result, except that the 
variation of shock strength with direction is not the sin 4 variation which Warren 
is forced to get as a result of his choosing to represent the lift as a distribution of 
doublets along the axis. 


If the aircraft or missile is accelerating along a straight path, the amplitude 
factor in equation (1) is modified and, further, if the path is not straight this modi- 
fication to the amplitude factor will also depend upon 6. Rao studied this and a 
brief account of his work, together with its application to a general winged body, 
is contained in Section 4. 


NOTATION 
r,9,x cylindrical polar co-ordinates 
x,y,z Cartesian co-ordinates 
U,, main stream velocity 
M,, Mach number of the main stream 


M local Mach number. In Section 4, M represents the flight 
Mach number of the body 


Mach angle 


U,.u_ perturbation to the main stream velocity in the main stream 
direction 


B= / (M,,?- 1) 
a,, sound speed in the main stream 


é local direction of flow 


y ratio of specific heats 
H,(x) Heaviside unit function 
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im R(x) radius distribution of the body of revolution 
- S(x) area distribution of the body of revolution 
S,,5.,5,.S, area distributions of certain “equivalent” bodies of 
revolution 
os d_ distance of the wing behind the body of revolution 
Z (x,y) thickness distribution of the wing 
on h (x)=Z (x, 0) 
he Z, (x, y)=aZ (x. y)/ax 
i a2 angle of incidence of the wing 
ck K,, K,,...K,... Bessel functions of the second kind with imaginary argument 
ng p Heaviside operator; p also represents the local pressure 
tal main stream pressure 
” Px density of the main stream 
ne (p-p..)/Px, pressure ratio 
en H(p) Heaviside operational representation of / (x) 
of U,.@ a perturbation potential 
U.o« the perturbation potential of a two-dimensional wing 
de U®, the perturbation potential of interference effects 
+ f.F functions used to describe the disturbance generated by the 
wing-body combination 
G_ function defined by equation (38) 
7(x,r) describes the characteristic curves 
T (x,r) describes the front shock wave 
(1(), 7) the parametric representation of a point on the supersonic 
leading edge of the wing 
x the angle between the plane of the wing and the plane 
through (/(»), 7) which is normal to the wave front 
ht s In Section 3, s measures distance from the wing in the plane 
through (/ (1), 1) which is normal to the wave front 
In Section 4, s represents the distance along a ray 
2,8 distance in the x-, y-directions measured from (/ (1), 
Ap/(4p.U.2) local lift coefficient 
k curvature of the path described by the body (Section 4) 
x’,y’,t variables of integration. In Section 4, t also represents time 
K,E complete elliptic functions of the first and second kind. 
a,b,c used to describe the shape of the wing plan forms in Section 5 
m=c/a 
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2. The Linear Theory 


As explained in Section 1, we start with the steady-flow problem; these results 
will be extended to non-uniform motion later. First the body and wing are treated 
separately, then these results are added, and finally the interference effects are 
included. 


Let the main stream velocity be U. in the x-direction, and let U..¢ (x, y, z) be 
the potential of small perturbations to the main stream. The function » must satisfy 
the partial differential equation 


» also satisfies appropriate boundary conditions on the body. 


For a body of revolution which has an area distribution S (x), and which is at 
zero incidence to the stream, the appropriate solution of (2) is 


x—Br 


0 


where the dash denotes differentiation with respect to t. This expression can 
be interpreted as giving the perturbation potential of a distribution of supersonic 
sources of strength U.S’ (t) which lie along the axis of the body of revolution. That 
the strength of the source distribution is U.S’ (t) can be deduced directly from the 
boundary condition at the body surface. This condition is (09/0r),-»=R’ (x), 
where R (x) is the body radius. 


Alternatively, the total change in flux due to the presence of the body must be 
U.oS, where 4S is the change in cross-sectional area. Hence U,6S/éx= 
2=U RR’ must be the source strength which satisfies the boundary condition at the 
body surface. 


The velocity potential for small perturbations due to a thin wing is the sum of 
the perturbation potential for a certain thickness problem and the perturbation 
potential for the problem of a flat plate at incidence to the stream. The thickness 
problem is symmetrical and the solution can be represented as being the potential 
due to a distribution of sources over the wing plan form; as for a body of revolution, 
the source strength is determined from the boundary conditions at the wing surface. 


For a symmetrical wing, whose plan form lies in the x, y-plane, the potential of 
the appropriate source distribution is 


AX, y, zZ= dx'dy’. (4) 


The source strength is 2U .Z, (x, y), where Z (x, y) is the thickness distribution of 
the wing and Z, (x, y)=0Z (x. y)/Ox. (x, y, z) is a right-handed system of Cartesian 
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co-ordinates. The double integration is over the part of the wing plan form which 
lies within the region (x x’)? > (y— y’)? + 


The flat plate at incidence is equivalent to a certain distribution of elementary 
horseshoe vortices over the wing plan form. The perturbation potential due to this 
vortex distribution is 


The strength of an elementary horseshoe vortex depends upon the local lift 
coefficient Ap/(4p,U..”). The range of integration is that part of the wing plan 
form which lies within the region (x— x’)? > B? (y- y’)? + B*z’. 


The perturbation potential for a thick lifting wing is the sum of expressions 
(4) and (5). 


If the combination of a body of revolution and a thick symmetrical wing at an 
angle of incidence is considered, the perturbation potential of the combination is the 
sum of the potentials due to the wing and body when separated, together with an 
interference potential. The interference potential corrects for the effect of either 
component at the boundary of the other. The interference effects of the body on 
the wing are O(R*h log R), where h is the wing thickness and R the body radius. 
This contribution to the interference potential is small and will be neglected. The 
effects of interference will be localised; they arise mainly from the interaction 
between the flow over that part of the wing close to the body and the flow over 
the body. 


In this region the wing is effectively two-dimensional. For a symmetrical 
wing the central section of the wing is given by z= — 2x +h (x), where z is the angle 
of incidence of the wing. The contribution of the incidence part to the potential 
can be shown to be at most O (2R*), which is small. Thus the main contribution to 
the interference potential arises from correcting for the effect of the flow over the 
wing at zero incidence at the surface of the body. 


In the upper half-plane the potential field for the two-dimensional wing is 
given by 


Ih 
oy = B (x- Bz). 


The component of velocity normal to the surface of the body due to this potential 
field is given by 
sin 0h’ (x— BR sin 4). 


The wings lie in the planes 6=0 and 6==. Here R=R (x) describes the surface 
of the body and h’(x)=dh(x)/dx. It is now necessary to find an interference 
potential ®; which will cancel this normal component of velocity on the body. 
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In cylindrical polar co-ordinates equation (2), which ®,; must satisfy, becomes 


or’ r or ax (6) 
This equation has the operational form 
1 oP, 1 
+ + 062 = B*p P,. (7) 
It is known that a solution of (7) is 
=1 


n= 


As the interference potential is symmetric about 46=47, only the terms for 
which n is even need be retained in (8). In equation (8) the A’s must be chosen 
so that 


= 


H (p) is the operational representation of h(x). From (9), 


= _-H(p) i -BpR sin 6 
4" | 
(10) 
—2H (p) i BpR sin 6 
and A,, BzK.,’ (BpR) sin cos 2né e ~ 3? dé 


BpR is small and, since K.,, (x) is O(x *") as x —> 0, it follows from (10) that A., 
is O(R*"*'h) for all n. 


Because all the K.,,’s have the same order of magnitude and the A.,’s are of 
decreasing order of magnitude at large distances, the most significant term in the 
expansion must be the first; thus 


2R (x) 


h’ (x) K, (Bpr)+O(R*hlog R). . ‘ ay 


When interpreted, the expression for <P; becomes 


(x, 


2 
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This interference potential is similar to the perturbation potential for a body of 
revolution. 


In combining the results which have been obtained, the axis of symmetry of the 
body of revolution is taken to be the x-axis, the y-axis to lie in the plane of the 
wing perpendicular to the x-axis, and the z-axis to be perpendicular to the x- and 
y-axes. The origin O of this system of axes is to be the intersection of the lead- 
ing edge of the wing with the axis of symmetry of the body of revolution. O is a 
distance d behind the nose of the axi-symmetric body. The body of revolution has 
radius R (x), where R (—d)=0; the area distribution of this body of revolution is 
(x)= 7R? (x). 


Then the potential field of the wing-body combination is given by 


x-—Br 
4n | V7 (x - x’) — B(y— - 


x-Br 
R(t) hr (t) It. a) 


0 
In equation (12), if —4R (t) h’ (t) is replaced by S,’ (#), S,(0) represents the area 
distribution of a body of revolution which, at large distances, has a potential field 
equivalent to the potential field of the interference effects. 


3. Correction of the Linear Theory 
3.1. BoDy OF REVOLUTION 


For a body of revolution which has an area distribution S$ (x), linear theory 
predicts 


u (14) 


x-—Br 
ay S” (1) 


= ox Vig 


0 


Here U.u is the perturbation to the main stream in the direction of the main 
stream, r is the distance from the axis of symmetry and x is the distance along the 
axis of symmetry. In this case the origin of the system of axes is the nose of the 
body of revolution. 


The detailed correction of the linear theory is contained in Whitham’s paper. 
However, the main points required in treating the flow far from a general wing- 
body combination are repeated here. 
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Correction of the linear theory requires that, in equation (14), x— Br be replaced 
by = (x,r), where s=constant is a more accurate representation of the characteristic 
curves, which in linear theory are the straight lines x— Br=constant. The new 
characteristics are curves because non-linear theory takes account of the effect of 
local fluid velocity and changes in pressure upon the speed with which disturbances 
propagate. In linear theory these effects are neglected and the speed with which 
disturbances are propagated is assumed to be constant. 


If now the result of replacing x— Br by = in (14) is approximated for r large 
compared with ;, we get 


— F (7) 
- om Ss” (t) 
where F (z)= | dt. 


0 


Further, because (p— p..)/p..= —yM,.°u, we get 


yM,,°F (7) 


Py 


Now along = (x,r)=constant, dx/dr=cot(u+<¢) where »=sin~' (1/M) is the 
local Mach angle and : is the local direction of flow. But 


y+1 
cot (u+2)=B+ - M.,? (v + Bu) +O (u? 
where v=09/0r; hence when dx/dr=cot(u+¢) is integrated, and the arbitrary 
constant chosen so that for small distances the corrected characteristic becomes the 
linear characteristic x - Br=constant, 7 is given by 


(y+ 1) (7) 


71/2 2 (17) 


7=x— Br+ 


3.2. SYMMETRICAL WING AT ZERO INCIDENCE 


The problem of a thick symmetric wing has been considered by Whitham”. 
At large distances we expect a solution of the form 


f (r, 8) will depend upon the shape and thickness of the wing in a way which is made 
apparent by a study of the full linear theory. 
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(42) 


FiGureE 1. FIGURE 2. 


In linear theory the perturbation potential for a symmetric wing which has a 
thickness distribution Z=Z (x, y) is given by 


The integration is over that part of the wing plan form in the region 
> B? (y—y’)? + 


The shaded portion of Fig. 1 is the range of integration of the double integral 
in equation (18). The curve C is part of the curve of intersection of the forward 
Mach cone from (x, y, z) with the plane z=0. 


The parametric equations for the supersonic leading edge of the wing are 
x=/(n), y=n. A point lies on the supersonic part of the wing leading edge if the 
tangent to the leading edge, at that point, lies outside the Mach cone from the point. 


In Fig. 2, AOCB is the leading edge of the wing, D and E lie in the xy-plane, 
i.e. the plane through A, O, C and B. CD lies in the main stream direction; CF is 
the intersection of the Mach cone through C and the wave front which is the 
envelope of all the Mach cones from points on that part of the wing leading edge 
which is supersonic. s measures the distance from the wing in the plane CFD, 
which is normal to the wave front and passes through the point (/(»),7). This 
plane makes an angle x=cos~' {/’ (y)/B} with the xy-plane. 


Whitham™ replaced x—/(n)- Bs by = and the Cartesian co-ordinates y and z 
by »+scosy and ssiny. A set of co-ordinates (x, 7, s) can be associated with 
every field point (x, y, z); » determines a point (/ (y), n) on that part of the wing plan 
form which is a supersonic leading edge. 
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C=(£(n)m) 


A FIGURE 3. 


If the origin of co-ordinates is moved to the point (/ (n), n), and the substitutions 
indicated are made, relation (18) becomes 


where z=x’—/(n) and B=y - ». 
If (19) is approximated for z/s small, it becomes 


In (20) the range of integration is that part of the wing plan form which lies in 
the region 2-/’8 <7; this region is shaded in Fig. 4. As 7/s—>0, 2-I’'B=7 is 
the limiting shape of the curve of intersection of the forward Mach cone from 
(x, y, z) and the plane z=0. 


For a wing with a subsonic leading edge the wave front is a Mach cone. The 
planes which are normal to the wave front all pass through the axis of this cone 
and the point (/ (7), ) is the vertex of this cone for all values of x, which, in this 
case, is equal to @. 


Now (20) can be written 


x) 
(2Bs) 
Here f(7.x)= dat 


The Aeronautical Quarterly 


whe 


Thi 


wri 


an 


an 


WING-BODY SHOCKS 


FIGURE 4. 


where the dash denotes differentiation with respect to ¢ and 


By 
Ay 


This integral is taken along the line x-/’8=tr; 8, and 8, are the ordinates of the 
points of intersection of the line 2 - /8=t and the boundary of the wing plan form. 


It follows that u= 


ax ~~ J QBS)" 


where a dash denotes differentiation with respect to >. The expression for u can be 
written 


=— F (=. x) 
J/(2Bs)’ (22) 
y= 1 
where F(r,x)= dt 
and the dashes denote differentiation with respect to 1. 
P-Dx, _ YM,,’F x) 
and 7, the corrected expression for the characteristics, is given by 
7=x—-—I(n)- Bs+ . (24) 
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From equations (23) and (24) and the form of the function F (-, x) it is plain that 
it is possible to find a body of revolution which at large distances produces the same 
disturbance as the wing. This equivalent body of revolution is centred on a line 
which passes through the point (/ (7), 1) and is parallel to the x-axis; it has an area 
distribution 


Be 
$,@=2 Z (2, B) dB. 
By 


The integral is that of (21) with Z (z, 8) as integrand instead of 0Z (2, 8)/0. 


For very large r. —> and s—>r. Hence becomes 


f (7, 


(2Br)* 


The expressions (22), (23) and (24) become 


FD) 
= ax 
P-DP. , F (7,9) 
Pee 


Here » is chosen so that /’ (y)=Bcos 4, and the equivalent body of revolution is 
assumed to be centred on the x-axis. 


The shock wave cuts off the characteristics before they intersect. If the shock 
wave is given by 7=T (r, #), the slope of the shock wave is given by 


dx (y+1)M, PG, 6) (r+ 1) 2 


where the dash denotes differentiation with respect to T. 


By introducing a further condition, namely that to the first order in the shock 
strength the shock wave bisects two intersecting characteristics, it can be shown that 


dx Mat F(T, 8) 
‘dr 25/2 Bs/2 pil2 
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Hence dr F’ (T, 4) r'! 


Multiplying both sides of this equation by 2F (7) and integrating yields 


(7, 9) dr 
(2B)? 
2 
Since F (T, #) is known from the shape of the wing, (26) and the last two expressions 
of (25) give the strength and position of the shock for large values of r. 


It is seen from (26) that, for very large r, T is near T,, where F (T,,9)=0. It 
follows that F (T, ) will be given approximately by 


Ty 


Thus, for large r, an asymptotic value of the shock strength is 


For distances of this order of magnitude the position of the shock relative to the 
linear characteristic x—1(n)- Br=T, is given by 


/2 
T, =x -1(y)- Br+ 21/2 ) 


Whenever the part of the boundary of the wing plan form which is a super- 
sonic leading edge contains a straight section within its length, the large-distance 
approximation (which involves approximating to a straight line the curve of 
intersection of the forward Mach cone from the field point and the wing plan form) 
does not yield an estimate of F (7,6) for small +. However, the range of 7 over 
which it is impossible to find F (7, #) in a simple way can be made arbitrarily small 
by confining the investigation to regions sufficiently far from the wing. In these 
regions 7 must be near T,, the first zero of F (r.#). Hence, for large distances, the 
values for the shock strength and its position which are obtained by applying the 
theory in the ordinary way will be all right provided that the straight section is not 
too long. 
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To find F (7, 4) when @ is given by I’ ()=Bcos 6=m, where m is constant, we 


integrate Z.=0Z (z, 8)/02 along lines which are parallel to the straight edge. This 


gives 


S, (0) and F(r,0)= - at. 


3.3. WING WITHOUT THICKNESS AT INCIDENCE 


A flat plate of the same shape and size as the wing plan form is inclined at a 
small angle to the main stream. At large distances from the wing the solution will 
take the form 


f(r, 
QBr)° 
In this case f (7,4) depends upon the distribution of lift over the wing plan form. 


The connection between f (7, ) and the lift distribution is made apparent by a study 
of the detailed linear theory. 


For a flat plate at incidence, linear theory gives 


l z(x- x’) 1 Ap 


(x, y. 


If the Cartesian co-ordinates are changed, as they were for the problem of the 
thick wing at zero incidence [(x=7+/(y)+Bs, y=y+scosy, z=ssiny, where 
Bcos x=’ (n)] relation (30) becomes 


J (cos x— 8)? +s" sin? x /[(Bs+7— 2)? B® (s cos — 8)’ — sin? x] 
Ap (2, 8) 3 
xT dadB. . (31) 


Again z= x’ -/(n), 8=y’—», and s and have the same meaning as before. 


If s is large compared with ;, (31) becomes 


| | V[r—(2-PB)] dzd8. . 


The double integration is over the part of the wing plan form which lies within 
the region z2—l’/B <r. 
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Equation (32) can be written in the form 


f(7.x)_ 
(2Bs) 
By 
oo 


The integral for S,’ (t) is taken along the line z—/’8=1; the ordinates of the inter- 
section of this line and the boundary of the wing plan form are f, and £,,. 


_ % __ F(x) 
we can write ax VQBs)’ (35) 
where F(7,\)= of 


P- Pw _ x) 


where F (7,\)= —G (7.x) J (37) 


As for the thickness problem, the improved equation for the characteristics is 


(y+ 1) M.* 


Bs+ Bale 


F (r, x) (39) 


Because of the form of G (r, x), it is evident that (36), (37) and (39) bear a close 
resemblance to similar relations for a certain body of revolution. The axis of this 
body of revolution is parallel to the x-axis, and it passes through the point (/ (n), ») 
where I’ (n)=Bcos x. This means that for a flat plate at incidence the perturbation 
potential at large distances is equivalent to the perturbation potential at large 
distances from a certain body of revolution. However, it should be noted that the 
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equivalent body depends upon direction; there is a different equivalent body for 
each direction. The area distribution of the equivalent body is S, (x) where 


S,0=S, (| F504), 


S, (x) is given by (34). 


For 0 <x <7, the perturbation potential of the inclined flat plate is minus 
one times that for the equivalent body. For - =<. <0 the perturbation potential 
of the plate is equal to that for the equivalent body. 


Now, again as for the thickness problem, if the shock wave is given by 
r=T (s, x), then at large distances from the wing the relation between T and s for 
any particular value of y is given by 


T 


F (7, x) dz 
For very large s, \ —> 6 and s —> r and (36). (39) and (40) become 
P- Px =" » F (7, A) 
1 
7=x-1(y)—Br+ F . (42) 
T 
j F (7, 9) dr 
3/2 
and _ 28) (43) 


~ F?(T,6) 


F (7, 4) is known from the lift distribution, so (41), (42) and (43) give the strength 
and position of the shock for large values of r. 


It is seen from (43) that, for very large values of r, T must be near T,, where 
F (T,. 9) is zero. From (43) it follows that F(T, 4) will be given approximately by 
Ty 


(2B)? f 
FT. 9= pit) FG Oar. 
0 


Thus for large r an asymptotic value for the shock strength is 


Ty 
0 
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The position of the shock relative to the linear characteristic x -/(n)—Br=T, 
is given by 


Ty 
1/2 
2B)" | F 0) dr 
In (45) and (46), 
. . . 7) 
_1f 


Now » is determined from I’ (n)=Bcos 6 and the equivalent bodies are assumed to 
be centred upon the x-axis. 


For —t <6 <0, substituting (47) and (48) in (45) and (46) gives 


| sin 6 | 2/2 


Mat | ar 
and T,=x-I[(n)- Br+ 


| sin 6 


(49) 


0 
Be 


When the supersonic leading edge has a straight section within its length, for 
which l (,)=m=B cos 9, then F (r, #) cannot be found for small +. However, as 
for the thick wing at zero incidence, at very large distances from the wing the range 
of s over which F (7, 4) cannot be found is small. Over this range F (r, 6) is taken 
to be the value obtained by applying the theory as if the edge were not straight; the 
error in doing this will be small if the distances involved are large enough. 


To find F(r,9) when Bcos6=m, we integrate Ap/(4pU,,”) along lines 
parallel to the straight edge; this gives S,’ (f) and 


room 


3.4. WinG-Bopy COMBINATION 


In linearised theory the expression for the perturbation potential of a winged 
body is given by equation (13) in Section 2. When this expression is corrected and 
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the large-distance approximation is made, it is found that 


dt+ 


-d 


+H, (r-1(n)) 


Bsind 1 f AO) 
2 


1@ 


dt (51) 


H, (t) is the Heaviside unit function and Il (n)=B cos 8. 
At large distances, if the shock is s=T (r, 6), 


P- Poo F (7, 4) 
> =yM..? 
V(2Br) 

T 
xe 9) dr 

F? (7,9) 


The relations numbered (52) locate the shock and give its strength. From the 
form of F (7,4) and these relations it is evident that in every direction the wing- 
body combination can be associated with some body of revolution. Further, the 
flow at large distances from this body of revolution will be the same as the flow at 
large distances from the wing-body combination. Because of its complex shape 
it is impossible to say very much about what happens near to a general wing-body 
combination. However, it is not necessary to use this theory up to the body for, in 
the regions close to the body, linear theory does not require correction. 


As r—> ©, it follows from (52) that asymptotic values for the shock strength 
and its position are given by 


1/ 
(2B)°!? j F (7, 9) dr 
P= _ 0 
ris 2B (y+ 1) 
Ty 
| ar 
(y+) Mx! 


x-1(n)-—Br=T, - pla 


4 
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4. Non-Uniform Motion 


When an asymmetric body accelerates along a curved path, further directional 
effects must be included. These directional effects are due to non-uniform motion, 
which alters the way in which the disturbance decays as it propagates outwards 
from the body. If the body is axi-symmetric, there are no initial directional effects; 
all the asymmetry is a result of the non-uniform motion. In this case Rao® found 
the excess pressure ratio to be 


P~ Po (7) 


— 1) 


where A= 1)'/?Mx« cos 6+ (a,,/ | dt) (55) 


In regions of compression, where F (7) >0, shock waves are formed; using 
body surface and are the orthogonal trajectories of the wave front. As for 
geometrical acoustics, the disturbance propagates along these rays. « is the 
curvature of the flight path. 4 is the polar angle in a plane perpendicular to 
the flight path, and a,dM/dt is the acceleration along the flight path. Ma,,’« is the 
component of acceleration perpendicular to the flight path. s(1-s/A) is propor- 
tional to the cross-sectional area of a bundle of neighbouring rays. ¢, is the time 
at which the ray set out from the body. 


Figure 5 illustrates the co-ordinate system used in this section. AB is the 
flight path, COD is a plane perpendicular to the flight path, where O is the position 
of the body at time t=f,. SCO isa plane inclined at an angle @ to the instantaneous 
plane of motion of the body at the time f=f,. s is measured along OS, which is a 
tray. OS is inclined at an angle 4x— , to OB, where »,=sin~'{1/M (t,)} and M (A) 
is the flight Mach number of the body. 


F (r) gives the initial “ wave” profile of the disturbance. Rao verified that, 
for an axi-symmetric body, the function F (7) in non-uniform motion is exactly that 
used to describe the flow pattern for the same body moving uniformly. 
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If a is the local sound speed and u the local fluid velocity along a ray, the 
correct speed of propagation of disturbance is 


*) =at+u 


In view of expression (54) for the excess pressure, this relation, when integrated, 
gives 


(+1) 


(t,) (t-—t, - s/a,.)+ {8 (M?-1))*? 


F(r)b(s), (56) 


s 


ds 
where b(s)= | al 
0 


When this integral is evaluated, we get 


b (s)=2A"/? sin-? (s/A)'/?, if A>O, 


and b(s)=2 (—A)'? sinh! {s/(-A)}?”, if A<O. 


In (56), the arbitrary function of + which occurs in the integration has been 
chosen so that near the body the improved expression for the characteristics 
becomes U (t,) (t- ts —s/a@,,). 


In regions of compression, where F (r) >0, shock waves are formed; using 
the Rankine-Hugoniot shock conditions, with the further condition that at the 
shock the disturbance propagates with the speed of the shock, it is found that, if the 
shock is given by 7=T (s), the relation between T and s is 


|F@a: 


M"? 
{32 (M? 1)} 


From (57) it follows that, if s is to be large, T must be near to T,, where 
F(T,)=0. Then, for large s, (54) and (56) become approximately 


Ty ) 
1/2 
J dr 
P Pew 2°/*M*!*y 0 
Po (M?- 1)'* (1 -/A) , 
2(y.+ 1)'/2M7!/4 1/2 
and T,=U (t,) 2 | Fe) ar} 
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The second expression of (58) locates the shock relative to the linear characteristic 
T,=U (t,) (t- t, —s/a,.). 


In view of Rao’s work, and the form of F (7,4) for a particular value of 9, it 
follows that the F (7,4) used to describe the shock pattern at large distances from 
the wing-body combination in the steady case can be used to describe the shock 
pattern in the unsteady case. 


Thus the position and strength of the front shock for the wing-body 
combination moving non-uniformly along a curved path are given by 


Ty 


[ Fo, 6) dz 
P— 


2 (y+ 
(M? - 1)'/4 


and T,=U (t,) (t- 


To 
F (r, 9) ar} 


F (7,4) is given by 


2x J J(r-2) 


(60) 


where H, (t) is the Heaviside unit function. 


In (60), 7 is chosen so that I’ (n)=Bcos @. (/(n), ») is the parametric representa- 
tion of the supersonic leading edge of the wing part of the wing-body combination. 
5. Examples 
5.1. A SYMMETRICAL WING 


In this example we find the variation with direction of the asymptotic value 
of the shock strength at large distances from a symmetrical wing, which has an 
elliptic plan form and a thickness distribution given by 


Z (x, y)= s(1- - 2). 


The thickness is zero on the boundary of the wing plan form. 
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From the theory, the disturbance at large distances is 


2F (7,9) 


By 
and \ = dp. 
By 


In the integral for S’ (t), z=x-I1(n) and B=y—»y. The integral is taken along the 
line z-/'B=t; 8, and £, are the intersections of this line and the boundary of the 
wing plan form. 4 is given by l’ (n)=Bcos 6, where, in this example, 


I(n)=a (1 = n? | 
85 (b? — 2abt 


and S” (t)= 7°) (b? +H ~s} 


88 (b? — n?)? 2ab 1/2 
yn { 
4/2 8b 


Hence F(r,@)=4¥2 [8 6,/a)?-51K (6,/a,)+ 


+18 - (61) 


K (k) and E(k) are the complete elliptic functions of the first and second kinds 
respectively. 


In equation (61), b,=7'/? and a,=(2ab)'!? (b? - 


(14 B*b? cos? 


Let = 


and let G(x)=(8x?-—5) K (x)+(8-—16x”) E(x). The first zero of G(x) occurs 
at x=0:56. 
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4/2 8b G (x) 


Th F (7, #= —— 
en (7, 4) 3n ( B°b? cos? aye 
1+ 


(62) 


Figure 6 shows the variation of G (x) with x. 
To find the position and strength of the front shock wave, we must evaluate 
T 
F (7, 9) dr, 


where 7", is the first zero of F (7,9). This gives 


x 
T, { xG (x) dx 


0 
_ 16/28b 


where X=0-56. 
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x 
| x G (x) was evaluated numerically to be 0:34. Hence (63) becomes 


Ty 
6) dr= 3/4 
0 (14 ) 


The asymptotic formula for the strength of the front shock is 


Ty 
P— Poo _ { (2B)'” 
= F (7, 9) dz . (64) 
Po r (+1) 
In this example 
P— Doo 0-97B'/481/2p1/2 


Poo cos? 6 3/8 
(14 


In Fig. 7, [(p- p)/P] (r/a)*'* is plotted against 6 for 5/a=0-087, b/a=2, and 
M=2. From Fig. 7, it is seen that, as r—> 00, at M=2, the shock strength in the 
plane of the wing is 0-358 times that in a plane perpendicular to the wing. 


In the direction 6=47, the value of shock strength predicted by the asymptotic 
formula for the shock strength is almost exact for r > 12a. 


188 The Aeronautical Quarterly 


of 
(25 


(F 
(p- 
equ 
equ 


is § 


axis 


mai 


May | 


0's 
0-4 
0°3 
4 
Py 
0-2 
O°! 
fe) 
90 70 50 40 30 20 10 (e) 
8 
FIGURE 7. 
0 
|_| 
|_| 
= 


WING-BODY SHOCKS 


When 6=0, the ratio [(p- p.,)/p.] to (p - p..)/ Po. is given for different values 
of r/a. (p-p.x)/ Po. 18 the exact shock strength and it is calculated from equations 
(25) and (26) by evaluating 


T 
F(T,0) and { F (7, 0) dz 
0 
(F (T, 0) is given by equation (61)) for a given value of T; r/a is found from (26) and 
(p—P»)/ Px, is the corresponding value of the shock strength given by the second 
equation of (25). [(p—p,..)/p.]’ represents the asymptotic value of the shock 


strength at the distance r/a, calculated from (26). [(p—p..)/p.]’ is given in 
equation (65). The resulting values are 


365-3 380°6 1460 9395 


(P Po)! Poo 
06298 0:7213 0-8281 0-9000 
[(P - Pool 


The position of the shock relative to the linear characteristic x—/(j)— Br=T, 
is given by 


Ty 


491/2 1/2,1/4 


In this example, 


(r/a)}!4 
Bb? cos? 
(i + 


2 1/2 
x= 1 Br=0-627a (1+ 


5.2. A WinG-Bopy COMBINATION 


Consider the effect of a wing-body combination at large distances from the 
axis in the directions 6= +47. 


The body of revolution is given by 


(8x/D (l- x), for O0<x<l 
0, for x<.0,x>1. 


The wing, which has a delta-shaped plan form, is inclined at an angle 2 to the 
main stream; it has a thickness distribution which is given by 


Z=8,(y) {1- } 


and =, (y)=6, (1 - y*/@’). 
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The wing plan form is illustrated in Fig. 8. The thickness of the wing is zero on 
the boundary of the wing plan form. In this example the vertex of the wing plan 
form is placed at a distance of one quarter of the body’s length from the nose. 


The contribution from the body of revolution to F (7, +47) for the com- 
bination is 


For 0<7,<0:75, the contribution to +47) from the interference 


potential is 


Here, and hereafter, +, =7 -0-25/ and c=0-75/._ The contribution from the inter- 
ference potential is zero for 7, < 0. 


The contributions to F (r, +47) from the thickness of the wing is found by first 
integrating 0Z (x, y)/dx along lines x=t, where ¢ is constant; this integration is 
between the points of intersection of the line x= and the leading edges of the wing. 
Next, we differentiate this integral with respect to t, divide by =/(r,—1f) and 
integrate the result with respect to t for O<:.t<7,. Thus, for 0<7,<c, the 
contribution to F (+,+47) from the wing at zero incidence is 


165, 


* 8 4+7,/0+ 


The contribution is zero for 7, < 0. 
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To find the contribution from the wing at incidence we replace 0Z (x, y)/dx by 
Ap/(4p,.U?) and proceed as for the thick wing at zero incidence; this result, 
multiplied by + B/4, is the contribution to F (r,+47) for the combination, due to 
incidence effects of the wing. 

For a delta wing swept-back inside the Mach cone, Goldstein and Ward have 
shown that 


Ap _ 4am,* 
p.U? E (k)/(m,? - y?/x?) 


E is the complete elliptic function of the second kind to modulus k= (1 - m,?B?) 
and m,=1/m. 


It follows that the conribution to F (7,+47) due to incidence effects is 


For 7, <0 this contribution is zero. This result is that obtained by Whitham™. 
For 7, > c this contribution to F (7, +47) decays in a way typical of a finite body. 


F (z,+47) for the combination is the sum of the contributions from the wing, 
the body, and the interference effects. 


)ror wine 
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Figures 10 and 11 are of F (7,47) and F (7, - 47) for the wing-body combina- 
tion, where the length of the body is /=1, c=0°75, 5,=0-01, 6=0:2, M=2, 
m=cot 29° and «=14°. Fig. 9 shows the variation with 7 of the individual con- 
tributions to F (r,+47). Under these conditions 


T) 


F (7, 4x)=0-0046, 


and | F (r,- 4x)=0-0161. 
These integrals were evaluated numerically. 


__ For a given Mach number the asymptotic value for the shock strength at large 
distances in the direction 6=4z is proportional to 


{freimar} 


T, is the least non-zero value of + for which F (r, 4x)=0. 


Similarly in the direction = - 47, the asymptotic value for the shock strength 
at large distances is proportional to 


{ [F (r, - dr} 


T, is the least non-zero value of + for which F (7, - 47)=0. 


For this example, when 6=4r, 


P-Px  0°084 


and when 6=.—4r, 


P—Px _ 0°156 


When M'=2 the asymptotic formulae overestimate the strength of the front shock 
wave, produced by the combination, by less than 10 per cent of its true value if 
r>91 when 6= - 42, and r> 2036 when 9= +42. The unit of length is taken to 
be the overall length of the combination. 
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Errata 
February 1958 (Vol. IX, Part 1) 
The paper “The Wave Drag of Wing-Quasi-Cylinder Combinations at Zero 


Incidence” by L. E. FRAENKEL contains the following printing errors :— 


(i) In equation (7), p. 59, t should be f on the right-hand side of the first line; 
in the second line, the limits 0,27 have been omitted from the first integral sign. 


> 


(ii) In Table I, p. 64, the ¢,,? in the top left-hand corner should be underlined. 
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